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SYNOPSIS 


M. AEUIIUOiM, Ph.D, 

Department of SlGctrxc''! Engineering 
Indian Institute of Tochnologji^, Kanpur 
October 1971 

OPTIMAL REGTJLATOES TOR SYNCHRONOUS MACHINiSS 

In recent times, tlic attention of power sysiem 
Engineers hce boon focussed on methods of increasing 
transient and dyntimic purlormance of power systems in 
order to minimize or elimiiiaue tho of foots of soverc 
system oscillations, usiiig modorn control theoretic 
concopts. By regulating tho terminal voltage and opeod 
uo some fixed roferenoo, it nas boon ohown tiuc the 
transient porf ormanco of the syncioro louo generators ocn 
bo improved* 

The conventional design procedure !'or suen 
voltage regulators and spood governors assumes apriori 
hnowledgo of suitable oonfigur'’tions for tho regulating 
and stabilizing equipments* This has drawbacks in the 
arbitrary choice of the controller configurations aiid 
the out and try procedure involTred in bhc c elect ion of 
regulator parameters to meet tno design sp»-cificatione 
such as overshoot, steady sbato orror and stability 
limit* Even when these specifications are met, ic may 
nob bo the best possible design. The notion of ruodom 



ann optimal contjTol bheory provides a convenient frame- 
work for incorporating bhe design specifications and 
for studying the system dynamic behaviour. Through the 
proper choice of an appropriate performance criterion, 
it becomes possible to impart the desired features to 
the system transient and dynamic behaviour. 

The ob]ecLive o' tlie present thesis is to 
formuHate bh' synchronous machine control as an optimal 
control problem and u'^'cn to obb-iin an integrated control 
for the excitation aiid piimo ^lovjr .‘^'''t-ol, '•<1 lOii Is 
optimal with respect Lo a chosen p^'^’rormance cj'ib r-'ion. 

Chapter I Is a gcaoral ixitroauc bion giving a 
brief discussion on bx^e synchronous machine conux-ol 
problem and bhe modern dtvelopment with particular 
reference to the present status of o.ttimal conbrol of 
power system dynamics. 

Chapter II is d'-'votod to the dovclopmonb of a 
mathematical model for Aiq aynoi^ronouo machixue co_'Uccted 
to an infinite bus system. The stabo space model oC "lie 
si33gle machine system is derived using winding currents 
as some of the state variables instead of j‘lux lixilcagco, 
in a form most suitable for the applioatiOwi of o 'C'-mal 
control theory. For small disturbances in the system 
which occur continuously during the normal opera uion, 
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overall configuration of the controlled systera. The 
transfer matrix relating the outputs to t^c inputs arc 
also derived . 

The observer designed in the last chapter introduces 
d 3 mamics in the control loops and also exponcmjially decays 
error. Henco it is desired to find a control law which is 
function of mcasurablo outputs. Chapter VI deals wit3', tti© ! 
suboptimal control of c nicfironous machine s^'stem. Hero Ihe 
regulator is oonsbr i ■« jd bo be a linear time invnx*iant 

I 

function of the mcasuriblo system outputs. Using an i 

itorativo algonthro uhc feedback control matrix is deter- 

1 

mined. The effect of ±lo suooptimal control law dctunninedi 

I 

for a particular opciating point when used at cliXfurwiit lo-j 

I 

levels IS discussed and conclusions are diavm, J 

I 

I 

The stability of a oon'^'entional synchronous maenmo 
delivering loading pov/er factor loads is poor beo. use of | 
lower excitation. Tho problem of maintaining stability unde 

these conditions by usirg a divided winding rotor sjmchronoi} 

1 

machine is discussed ii- Chapter VII, A st‘‘to space model is 
derived and the behaviour of the system with c cjiVcsAitional 1 

I 

angle and vdltago regulators is invcs uiga cod . I 

i 

In Chapter VIII, tho optimal and subopt imsl contiol [ 

1 

of tho divided winding rotor synchronous ivchino u dealt ^ 
with. The system bch.viour with the optimal cc .trollor ! 



for largo disturbances is discussed. The superiority 
of optimal regulators over conventional regulators is 
osLablishcd . 

Tho disturbances occuring in a power system are 
random in nature and also uncorbainties are introduced 
in the output measurements. Hence the design of optimal 
regulators should take in bo effect tho presence of 
random disturbances. In Chapter U. the opbimal control 
of synchronous machine in uho presence of such renadoin 
disturbances is described. The aisturbances are assumed 
Lo bo wliito noise with knovm soutistics. The average 
aciiriviour of the syotom in thv. prcKjeiice of noise is then 
dcbermined. The optimal control law for those conditions 
IS bhc same as that for the doterminis bic case, but the 
porformanoo index value is increased on an average* 

In response uo the demands of tho system gro\7th 
and more interactions of control systems, chc dynamic 
analysis becoraos more expensive and time consuming oven 
on a fast digital computer. Tho question mitjht, bhereforo, 
be asked whether simplified models might not bo app2?opriabe 
for suoh systems. A method of simplify’‘ing such large 
linear dynamic systems using Schwarz canonical form which 
does not require tho computation of cigon values and 
cigon vectors as hithor to done in existing mebhods, lo 
described in Chapter Z, Tho response of the original 
and simplified models is compared. 



xix 


The nctiiod of stotc variiblo grouping is used to 
obt_in c3 Ti^Ouccd mode] tor nulti variable control systems. 
Using tbe IV dueled nodclj an optinal regulator is obtained. 
Tills Cl ntvol lav' is used as a suboptinal control lav/ for tlic 
original system. The pex-formanco of the original system is 
investigated with this suboptimal control law. 

In the concluding chiijter, the results of the 
thesis are reviewed and future lines of resoaren are 
delinatcd . 



CHAPTEPl I 


nriRODUCTioF 

1 .1 SYITCHEOHOUS MACHIEB COlTTEOIi 

In the developing countries of the world, the use 
of electricity h^s been doubling every 10 years, Ihe 
rapid and continuing expansion has led to the large extra 
high voltage interconnected power systems spreading 
across national boundaries. Ihe growing demand for 
becter security and contiAiuity of service on one hand 
and investment and operating cost reduction on the other 
hand affect the role of controls both in system de.^ign 
and operation. Even though the security is uxie most 
important aspect of q.uality of service, the ccaitrol 
systems play an important role in minimizing the syscem 
frequency and voltage deviations. 

With advent of fast acting voltage regulatiors 
and significant development in control system theory, 
the study of stability of synchronous machine sys-em 
provided \^ith these regulators, calDed the dynamic 
stability, IS assuming much importance, fhe syncimror.ous 
machines provided vath continuously acoing t,utom«-uic 
voltage regulators can be operated beyond tneir steady 
state stability limits, Ihe machine falls out of step 
or becomes unstable wnen the dynamic stability Ixiiit is 
reached f Ihe limt depends not only on the machine 



QjifiLPSdJsyi'S'bios biiij 8lso upon "bli© systsoiu losd snd 

characfeensliics • G-reater reliance iSj1iherefore» 
placed on well designed control systems. Automatic control 
tecliniqu©® Q^a conceived as a means to obtain better 
operating performance for a given system and these are 
used towards improvement of overall system response by 
designing proper feedback loops. The concepts of modem 
control theory and methods from multivariable system 
matrix analysis and optimization have a very important 
bearing on the basic engineering approach to the system 
analysis for design and operation. 

In the conventional methods of design and 
analysis of synchronous machine for control of voltage 
anfi freCj.uencyi the machines were represented by simple 
models* The oon1a?ol of voltage and freq.uency is viewed 
as two distinct control systems • Separate configurations 
are chosen with no interaction between them* The freg,uency 
is controlled by prnaemover governor which derives a 
feedback signal from the rotor speed and its derivative . 
More recently rotor angle control is also used. Similarly > 
the voltage regulator derives a signal proportional to 
the terminal voltage and/or some function of the same. 
Apriori configurations of the regulators for both 
frequency ani voltage controls are assumed independently. 



The complete model including the machine dynamics 
and regulators are derived as a single higher order differ- 
eniial equation* The time constants of the various 
rc gala ling and stabilizing equipments are assumed and then 
the gains of Ihe amplifiers in the main and stabilizing 
loops are determined from stability considerations. The 
transient performance specifications cannot be incorpora- 
ted directly into this design. Once the gains and time 
cons bants are solccted, the system is analysed for 
design specifications such as settling time, steady state 
error, TnayTTrum overshoot etc. The dooign is a trial and 
error process because the regulator parameters are 
adjusted till the design spcoi.”ications are nearly mot 
with , 

The methods used for bhe analysis and design of 
classical control systems are the Routh-Hurwitz and 
Hyquist criteria and root locus ■cochnique , These methods 
are frequency domain analysis oriented. The specifications 
must be given in terms of frequency domain oharacterisbics. 
More recently, the mebhod of sensitivity analysis and 
Microvio method are also employed to obtain botbor 
regulacor parameters. The classical control techniques 
can be easily applied only for small systems and tliat too 
for linear and time invariant systems, These methods cannot 
be easily extended to multivariable control systems. One 
cannot obtain the best poosible design by these methods, 
since controller oonXiguro lions and the feedback signals 
|are fixed apriori. 



1 ,2 ISaDERir TEEBDS 

With the growth of system size, the analysis and 
design of control systems become more difficult without 
the use of digital computers* A new approach to the 
synchronous machine controller design using optimal 
control theoretic concepts is increasin^y used* The 
synchronous machine control is identified as an optimal 
control problem* An integrated form of control law is 
obtained for both frequency and voltage controls. Since 
the system behaviour is governed by all the system 
abates, it is necessary to take into account the changes 
in all the variables and thus derive control signals 
which are functions of all these quantities. The control 
of voltage and frequency is treated as a siaagle control 
problem. 

In con^jarison with more conventional methods for 
feedback control system design, the modem optimal 
control theory design procedure is more direct because 
of the inclusion of all the mgjortant aspects of perfor- 
mance in a sin^e design index. Also the method can be 
used to nonlinear as well as time varying systems, however, 
at the e 23 )ense of increased coi^putational complexity. The 
conversion of prescribed design specifications into 
meaningful mathematical performance index may not be 
straight forward. But a number of performance indices 
are suggested in the literature and the nonuniqueness 



of bhem makes the selection of a performance index 
simpler. The modern control methods require complex 
computer programmes and a good deal of computing time 
for nonlinear and time varying systems. This is not 
a disadvantage in itself because the classical control 
methods cannot be used as such for the nonlinear and 
time varying systems. 

An optimal regulator obtained in reference 18 
for tne single machine system uses a simple model for 
Lho machine and flux linkages as soone of the stace 
variables. The configurations of the governor and 
voltage regulator are fixed apriori and the inputs to 
those regulating equipments are obtained by optimal 
control theory. Also it does not investigate the 
suitability of the optimal control for large disturbances 
and it is assumed that all the state variables can be 
measured for feedback. The reference 20 considerg a 
current model for the machine and the optimal regulator 
performance is investigated for step load changes. The 
effect of incomplete state feedback was analysed for 
impulse disturbances in reference 29, 

The progress in power semiconductors has made 
static excitation of synchronous machine possible using 
controlled bridge rectifiers. By fast electronic control 
all time lags except the mam generator field time 
constant are praotically eliminated « Blectrohydraulio 



governors combined with electronic controllers provide 
better frequency control and improved dynamic performance 
with lesser governor dead bands compared to conventional 
mechanical or hydraulic regulators. Ihorefore the 
dynamics of -the voltage regulator and speed governor 
need not be considered along with the machine model. 
This thesis makes a detailed study of the design of 
optimal and suboptimal controllers and investigates the 
performance both for small and large perturbations* The 
physical realizability of the controllers is also 
discussod. 

1 ,3 SOOEB OP THE THESIS 

The synchronous macliine control is identified as 

an optimal control problem. I’or successfully applying 

the techniques of optimal control theory, a state space 

model is derived in a suitable form for the synchronous 

machine system. The resulting optimal control policies 

for nonlinear systems are difficult to implement. Hence 

a linearized state space model is obtained for which a 

21 

constant feedback optimal control law is derived . A 
fairly common type of voltage regulator and a speed 
governor are chosen and the performance of the system 
is investigated with these regulators. It is shown that 
the performance is affected largely by the regulator 
parameters. An improper choice of them may lead to 
system instability. The conventional methods of 



selecting such parameters are discvissed* The second 

mechod of Lyapunov is employed to determine the optimum 

1 6 

values of the regulator gains y using a state space 
model for the machine and regulators. 

Prom investigation of the system performance 

provided with optimal regulators for large disturbances, 

it IS concluded that the optimal linear feedback control 

law can be used for tho nonlixiear system. The implemen- 

tvition of such a control law requires the direct measure' 

mont of complete state vector. It is seldom that all of 

Inem can be monitored directly and easily. To overcome 

chis difficulty, a compatible dynamic obsoirvor is 

designed vdiich will rccon struct the state vector from 

the measurable outputs . However, it is found bhat the 

cascaded optimal control system with the observer can be 

used only for small perti rbationc because the overall 

feedback control law is dif Joront for different 

operating conditions. Therefore, the question of using 

a control law which is a linear feedback of the output 

28 

ViU'iables IS discussed . The suboptiraal control law 
thus obtained can be used even for largo disturbances 
and hence it can be easily implemented on bhc nonlinear 
model. 

The normal synchronous machines have limited 
reactive power capacity and. the stability of the system 
under lightly loaded conditions is poor. The reactive 



capacity cannot be unproved even by fast acting voltage 
regulators* By the use of an additional field mnding 
with proper control schemes, the reactive limit can be 
incrocsod as well as the dynamic performance of the 
system. A divided winding rotor synchronous machine is 
discussed which will improve the system performance. An 
optimal state regulator is obtained for the linearized 

34 

model and as such it can be used on the nonlinear system. 
The design of an observer and incomplete state feedback 
of the d.w.r. machine are invesbigated . It is found that 
the suboptimal control law gives a satisfactory perfor- 
mance even for large disturbaiices and therefore it con be 
practically implemented. 

The stochastic optimal control is an active field 
of research in the recent times* The synenronous machine 
control IS investigated in the presence of random distur- 
baiioes in the state and uncertainty in the output measure- 
ments. The effect of noise in the measurement of output 
variables on the overall performance of the system with 
opTiimal control is also studied. 

For multivariable control system design and 
analysis by the modem optimization techniques, excessive 
computer memory and tune are required. A simplified model 
IS obtained using Schv/arz canonical form for the analysis 
of large I linear and time invariant systems^^ . The 



fea&±bility of using reduced models to obtain suboptimal 
conui oilers for the complebo system is also investigated. 

1 .4 Ani OP IHE THESIS 

The optimal control bhcoretio concepts are 
iiicre 'siaigly applied to the power system design and 
cOiibrol in the recent years* The need for the study 
of synchronous machine voltage and frequency control by 
optimal regulators becomes important in the light of 
system complexity both in design and operation, The 
aim 0 ? tho thesis is to obtain a mathematical modaJ 
fOx olio single machine ays Lem in a suibable form nad 
then obtain an opLiiiWl rcgulacor ’'Inch can be easily 
implemented on pr-actic-l s-suons. large interconnected 
poi/er systems demaiid complex co iputation and iiiereforc 
the question of using simplified models for •'nalysis 
becomes rolcvtint. The opbimal control lay obtained 
for tho reduced system v/111 bo suboptimal for tho 
original system. Even if the control law is suboptimal 
the implementation of the same on pracT:ical systems 
should be feasible. These aspoebs are discussed in 
the coming chapters. 



cHiCPm II 

STATE SPACE MODEIi OF SYNCHEOEOUS MACHINE 

2.1 lETEODUCTIOE 

In this chapter, the state space model of a single 
machine system is derived in a form most suitable for the 
design of optimal regulators. The study of dynamic 
behaviour of synchronous machines has followed tradition- 
ally the analysis of the piecewise linear model using the 
characteristic equation of the system by employing classi- 
cal control techniques. A simple alternative in keeping 
with the modem system analysis is to formulate the system 
behaviour, not in terms of a single highworder differential 
equation but in terms of sets of first order equabions. 

Such a mathematical model is most suited for the application 
of modem control theoretic concepts and system optimization. 

A conventional voltage regulator and speed governor 
configuration is chosen and the state model of the controlled 
system is then obtained* The performance of the free and 
controlled system is obtained for an impulse type disturbance 
on the system. The responses of the free and c anventionally 
controlled systems are compared. 
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2 ,2 G-ESaEIUIr I019IiINEAIl SISTM £^UAO}ICMS 

The sxn^e machine system shown in Figure 2,1 

consists of a salient pole synchronous generator connected 

to an infinite bus bhrough a long transmission line. The 

state space model of this system is desired. The inductances 

of the various machine windings are functions of the 

instantaneous angular position of the rotor and hence the 

differential equations describing the machine performance 

2 

are time varying equations. Park’s transformation is 
applied to the stator quantities to obtain time invariant 
differential equations. The assumptions made in the 
development of the syscem equations are : 

1 . Saturation and hysteresis in every magnetic circuit 
and eddy currents in the armature iron are neglected. 

2. The stator windings are sinusoidally distributed around 
the air gap so far as mutual effects between them are 
concerned. 

3. Blectrical transients on the transmission system are 
neglected and also the line parameters are ass-umed to 
be invariant with instantaneous system frequency. 

4. Mutual inductances between armature and rotor circuits 
on the same axes are made equal by the proper choice 
of rotor base quantities, 

5. Only two damper bar winding are considered, one in 
the direct axis and the other in the quadrature axis. 



6. OjoCLy those modes of operation which do not re(iuire 
zero axis variables are considered. 

The geometrical configuration of different windings 

IS shown in Pigure 2,2. The performance equations are 

3 4 » 5 6 

given below for three phase balanced operation^ • All 
the quantities are m p.u. except time and angle which are 
in seconds and radians respectively* The per unit angular 
frequency is chosen as unity and hence reactances are used 
in the flux linkage equations, instead of inductances* 

Direct Axis Dlux Llnlcages: 

’*'fd ® ^ffd ^fd “* ^ad ^d ^ad ^kd 

’^d * ^ad ^fd - *d =^d + ^ad ^kd (2*2) 

*kd ® ^ad ^d •" *ad ^d ^kkd ^kd (2.3) 


Quadrature Axis Dlux linkages: 

’"4 = - *4 =^4 ^*4 ^i :4 

^4 “ ” ^4 ^Kfcgt (2.5) 

Direct Axis Voltages: 

*^f d “ Wq P '^fd ^fd % (2*®) 

^d *^P^d -^a^d (2-7) 

0 « Jr p ^kd + ^kd ^d (2‘®^ 



Generator 




1_4 j 

^ Vo ,5 

□— J 

— □ — 



Transmission line 

Infinite bus i 


FIG. 2-1 SINGLE MACHINE SYSTEM 



.m/ ...LJ3 


-/tfW-- 


FiG.2-2 SCHEMATIC DIAGRAM OF SYN. GENERATOR 





'te -liii' i 'w A 'M ' . -Of S, Y .N » /G E N £ R AT OR 




Quadrature Axxs Voltages: 


^ ^ *4 - ^ ♦d 

Electrical Torque in the Air G-ap: 

®e =* ^q “ ’^q ^d 
G-enerator Terminal Voltage: 


(2.9) 

( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


The bransmission line equations are obtained referring to 
the phasor dia;,ram (ri:,ure 2,3) as 


■V/j s= V^ Sin 6 + r. 1 j - i„ 
do e d e q 

V ss V Cos 5 + X ij + r 1 
q 0 ® ^ '“■e d ^ e q 


(2.13) 

(2.14) 


Dynamxo Squation of Mo1;ioxl of Hobor: 

a « 


“ $ = \ - '^a a^- 


(2.15) 


The variables used in the above equations are explained in 
the nomenclature. 


The set of equations (2,1) to (2.15) are manipulated 
to obtain the state spaoe model. The state variables are 
Chosen as the winding currents in the d and q axis of the 
stauor and rotor, the rotor angular velocity and rotor 
angle. The winding currents are ohosen as some of the 
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state variables instead of winding flinc linkages because 
this seems to be a natural choice when machines are to be 
interconnected through a transmission network where the 
performance may bo more readily visualized m terms of 
voltages and currents ra uher than in terms of flux linkages. 
The implementation of optimal regulators requires the 
measurement of state vector for feedback. Thus the current 
variables' offer definioe advantage over flux linl>.age 
variables as the former variables aro comparatively easy to 
measure than the later ones. Tho machine terminal voltage, 
rotor an^e and rotor angular velocity are assumed lo be 
the measurable output variables. The concrol variaolos arc 
selected as tJio input torque T^ and tne fiel(5 voliugo 
referred to the stator side. Substituting cur-XiiCs for 
the flux linkages from equations (2*1) to (2.5), the 
following equations are obtained : 

p 6 = w - Wq (2. 1C) 

pw =[T^ - K^(w - w^) ~ x^^ if^ 1 ^ - x^^ ij,^ 1 ^ 


+ (x^ - ^q^^d ^q ■’ ^aq ^kq 


(2.17) 


^ffd P ^fd ~ ^ad P ^d ^ad P ^kd * x^'J ^fd “ ''^o ^fd 

... (2.18) 

"^ad P ''fd " *d P ^d + ^ad P ^kd = ''a =^d ''o ^o ^ 

+ ’^0 ^e ^d - % ^e ^q + '^^Sq ^kq “ -q ^q) 

(2.19) 


« • « 



^ad " *aa P ^a ^Mcd P ^ka ■ " *0 

- p P ^k4 = "0(^0 "^a + ^e ^1 

+ r^ 1 ^) - w(2^ - *d ""a + ^a \a> 

... (2.21) 


- ^au P ""a + ^kfci P "hci = - "0 ’^kq ^4 (2.22) 

« = « (2.23) 

W = W (2.24) 

4 = (T„ Sin « + ij - 1^)2 + (Y^ Ooa « 

+ *e ^a + (2.25) 

The nonlinear s cato nodcl of tno ej’^stem cai> be obtained 
from equations (2.16) to (2,25) us 

Bis f(X) + D u (2.26) 

Y s s(X) (2.27) 


wliere the state vector X, tno ouupub vector Y and ulio 
control vector u are given by 

Z = ( fi , w, 1^^, 1^, ij^^, ig^, 

Y s ( « , w, V^)® 

tl s= ®1^^ 


The matrices B and D are giiren by 





0 

0 

0 

0 

0 

0 ' 


0 

1 

0 

0 

0 

0 

0 


0 

0 

^fd 

-^ad 

^ad 

0 

0 

Sr ss 

0 

0 

^ad 

-^d 

^ad 

0 

0 


0 

0 

^ad 

-^ad 

^kkd 

0 

0 


0 

0 

0 

0 

0 


^aq 


0 

0 

0 

0 

0 

-=^aq 

^kkq 


*" 0 


0 Wq 

W^ad 0 

0 

0 6" 

D B 









0 

^As■ 

0 

0 

0 

0 0 


Ihe Teotors f(Z) and g(I) are given by 


f(I) 


w - 

r- K^(w - Wq) - + 

^ag. ’*' ^^d " ^q} ^d ~ 

- “^o ^fd d 

Wq Vq Sin 6 + w^(r^ + r^) x^ - (w x^H^ 

* ^ ^aq ^q 


^0 ^kd ^d 

^0 ^0 « + (Wq Xg - w x^) + w^ir^ + r^)!^ 

- =^ad d - ==ad ^d 

*■ ''’^o ^kq ^kq _ 



g(x) 


6 

w 


[ CV„ Sin 6+ Tg 1 ^ - x„)^ + (V^ Cos 6 


e 


I 

L. 


+ 1^ + r^i„) 3 

e a e 


2,3 IiINBimiZBD STATE SPADE MODEL 

Except for a very few special cases, there are 
no exact analytical methods for analysang nonlinear systems. 
Practical ways of solving nonlinear problems involve either 
graphical or experimenbal approaches. Another method of 
analysis is the piecewise linear approach. The analysis of 
linearized model has been extensively studied in the 
literature, Purthor the rcsulbs obtained from the study 
of linear model can be extended with some approximations 
to the nonlinear case* It is desired to obtain a linear 
feedback control law which is easy to implement either on 
the linear model or nonlinear model. To apply the 
techniques of opuimal control theory it is convenient to 
represent the system as a linear time invciriant system. 
Hence in this thesis, the study of linear system is 
considered. 

Using Taylor senes expansion of equations (2.26) 
and (2,27) about an operating point and retaining the first 
order terms in the expansion, the state space equations 
reduce to 
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EZ = EX + Du (2,28) 

Y = 0 X (2.29) 

y/here X, Y and u are tibo deviations of the state, output 
and control vectors respectively about the given operating 
point considered . The matrices E and C are given by 



0 


1 

0 


0 


0 

0 

0 


0 

• 

-%Ai 

"^ad^ 

qoA 

^24 


'ad^qoA 

^26 

“^aq^dc/ 


0 


0 

-Wor 

fd 

0 


0 

0 

0 

r » 

W V 1 

”o 0 

C!os6q 

CVJ 

0 


0 


%^aq 

' 

] 

0 


0 

0 


0 


“•^o'^kd 

0 

0 

1 

-V 

oS^'o 

CM 

-Vad 

*o(^d+*e^ 

"^^ad 


0 


_ 0 


0 

0 


0 


0 

0 

-Vka 


” 1 

0 

0 

0 

0 

0 

0*^ 




C s: 

0 

1 

0 

0 

0 

0 

0 





J31 

0 

0 

°34 

0 

°56 

0 





where 

^24 * ^ ^^d ” ^qo ^aq ^qo ^ ^ 

*26 = - l^ad ^fdo + ^ad ^do + ^^d " ^g^^do > ^ 

*42 “ ^aq ^kqo ” *q ^qo 

*62 “ *d ^do ■■ *ad^^fdo ■*■ ^kdo^ 

®31 “ ■^o^’^do °°® 'o - ^qo ^o^Ano 
®34 - (’do + ’qo ’'e^Ano 

®J6 ' (’qo - ’do ^‘e'Aao 


The equation (2,28) is premult iplied by the inverse of E 
to obtain the system state node! as 

i = A I + B u (2.50) 

where 

A = F and B = E’^ B 

Thus the state space model m the required form for the 
single machine infinite bus system is given by equations 
(2.29) and (2.30) which are used in the coming chapters. 

The system chosen has the following parameters 
where all the quantities are given in p.u. except time ard 
angle which are in seconds and radians respectively. 

Synchronous Machine Parameters^; 


*08 = ■' -O 

^aq * 

0.6 

Xd « 

1 .2 

^q 

« 0.8 

f d* ^ 

^lclcd== 

1.1 

^Idcq® 

0.8 

^al 

= 0.2 

II 

o 

• 

o 

II 

0,02 


0.04 

^fd 

=: 0.0011 

Wq » 314.0 

M s 

0.0192 


0.0032 

^0 

= 1.0 


Transmission line Data: 

Tq » 0.05 Xg ss 0,3 for both lines 

The operating point data, when the machine is delivering 
rated iCVA at 0,8 p.f. lagging to infinite bus, referring 
to phasor disgram of Figure 2.3, are obtained as 

Sq * 26.3° s 0.8921 i^j^ s: 0.452 i^^^ » 2,2614 

,0 = 1-^864 = 1 .2377 W ^20 


0.3526 



For the above operating conditions^ the matrices A, B aM 0 


are calculated and are given by 


A s 


0.00 noo 0,00 
0,00 - 0.1667 '- 25.54 
- 541.18 0,6953 - 2.13 
- 1136.48 1 .4601 - 0.66 


0.00 0.00 0.00 0.00 
9.42 - 23.54 - 99.20 27.88 

- 36.23 24.15 664.23 - 362,31 

- 76.08 - 12.08 1394.88 - 760.85 
- 56.23 - 38,65 664.23 - 362,31 


-541 .18 * 0.6953 1 .33 

397.98 3.4027 897.14 - 1345.71 897.14 - 53.83 - 26.91 
298,48 2.5520 672.86 - 1009.29 672.86 - 40,37 - 35.88 


T 



“ o.oo 

0.00 

2.1255 

0,6642 

- 1 ,3285 

0,00 

0,00 

B ss « 

0.00 

52.08 

0,00 

0,00 

0.00 

0.00 

0 . 0 £ 










“ 1.00 

0.00 

0,00 

0.00 

C.OO 

0.00 

0 . 0<T 

G » 

0,00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 


- 0,17 

0,00 

0,00 

0.32 

0.00 

- 0.04 

0 . 00 _ 


By checking the real part of the eigen values of the system 
matrix A, it is found that the free system is stable. The 


response of the free system obtained by digital simulation, 
for an perturbation in the field current of 0,05 p.u. 

is shown in Figure 2.4. The response of the system is rather 
poor because it is too oscillatory with more overshoot and 
takes more than 10 seconds to settle down. The response can 
be improved by having proper controls. 



2,4 PERPOEMMOB WIOH OOIWEHE lOBAI. EEGULAIOES 




In the conventional design procedure of regiilsctors, 
suitable configurations for the regulating and stabilizing 
equipments are selected apnori, Then the system is 
analyzed assuming various time constants, to determine the 
overall loop gain from stability considerations. This 
gain IS distributed to various components in the control 
loop depending upon their physical feasibility. Ihe gains 
and time constants are further adjusted until the system 
specific at 10318 suoh as overshoot, settling time etc. are 
satisfied « 

A fairly common typo of vdtage regulator and speed 
governor are shown in Pigures 2,5 and 2.6 respectively. 

She voltage regulator has a stabilizer with gam Kg ai3d 
time constant Ig, for the amplifier m tho main loop. The 
speed governor has two time constants Tg and Tj^ representiiqg 
the governor and turbine time constant. The performance 
equations for the regulators are given by^ 


Voltage Regulator: 

p %a “ \ ) - T; % (2.31) 

P^B (2.33) 


Speed Governor* 
K. 

^'i 


2 m 

P 


— g 

"o g'^h g h ^ 


1 


T. 


■g*h 


* • s 


(2.33) 








The ato-VB equatioiis (2,31) to (2.33) describing the operation 
of controlling equipments can be linearized and added to the 
S3’'3tem equation (2.30). Choosing the state variables Xg, Xg, 
X^Q and as AVg, aT^ and I»aTj|. respectively, the 

state equations for the regulators are given by 


Xg — a-j(oj^ X^ + c^^ X^ + o^g Xg + Xg) - ^ X^ (2,34) 
4g =s + O 36 ^6 ^ ^ 8 ^ 



+ (aj - i^)Zg 

( 2 . 35 ) 

iio = 

^11 

(2.36) 


K 


iii . 

®3 ^10 

(2.37) 


whore 

^ 1 
^ ag « ai ^ a3 « - 5^ 

Combining controller equations (2,34) to (2,37) with system 
equation (2,30) the state space model for the conventionally 
controlled ^stem becomes 

^ as X (2,38) 

C 

The following regulator paramotors are chosen^ as 


» 5,0 

Kg * 0,04 

K . 20.0 


If 

« 

0 

Tg 0.5 

Tg » 0.5 

= 0.5 


The system matrix A-j for the operating point considered in 
the last section »ath the above paramotors is given by 
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0?he controlled system is again solved for the 
performanoe when there is an iir£^<4<L.£ disturbance in the 
field current of 0.05 p.u, The response obtained by 
digital simulation using Eunge-Kutta fourth order inte- 
gration method is shown in ligure 2.7. The response is 
still oscillatoiy and takes about 8 seconds to settle 
down. An improper choice of the regulator parameters 
may sometimes render the system unstable. 

2.5 OOJTOIiUSION 

An accurate stabs space model of the synchronous 
machine system is derived \Thich is most suitable for the 
application of optimal control methods. The systematic 
method of obtaining the state space model is indicated. 

The configurations for the voltage regulator and speed 
governor arc selected and the porformanoe of the regulated 
system is compared with unregulated system. Since the 
regulator parameters hove marked effect on the system 
performance, proper parameter values are necessary. 
Methods of selecting those parameters are discussed in 
the next chapter. 



CHAPTSE III 

OPiniAL ilEGDLAPOE GAHIS BT SBCOICD 
ME5THOD OP IiYiLPUlTOV 


3.1 IITIEODUCIIOH 

Ihe optimum selection of voltage regulator and 
speed governor gains* is discussed by using Iiyapunov*s 
second raccbod so that the set cling time is ininimum m the 
event of system disturbance. The voltage regulator and 
speed governor configurations are chosen aprion. 

Diffox’ont methods of analysis and design of synchronous 
machine regulators using conventional control techniques 
are discussed. The performance of the aystom with 
optimal regulator parameters is then obtained for 
comparison with the response given in Chapter II. 

3.2 raTIIODS OP DBSIGII UD AEAI.YSIR OP COJTVEHTIOITAL 
EBGllIiATOHS 

\/ith increasing system voltages and generator unit 
size for economical reasons, the loss of stability in the 
dynamic region of operation of any one machine becomes a 
more serious problem. Properly designed excitation systems 
with voltage regulators of the continuously acting type can 
raise the stability limits considerably. Greater reliance 
IS therefore placed on well designed control systems which 
extend the operating margin of staoility. 
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Once the mathematical model of the system including 
the alternator and regulators are established, the classical 
control methods can be applied for the analysis of the 
behaviour of the system with regulators. Since the develop- 
ment of 2-a3cis theory by R.H. Park in 1929, the general 
effecu of regulators on the stability limit under particular 
operating conditions have been discussed by many authors. 

The nolliods used for the selection of regulator parameters 
which give a better performance are Routh-Hurwitz and 
ITyquist criterion, root locus techniques, sensitivity and 
Mitrovic methods. 

Concordia used the generalized machine tjoieory for 
modelling Lho system and the well known Routh’s criterion 
to dotermino the stability limit of round rotor machines 
under unity power factor operating conditions. 

8 

Messerly extended the analysis to include both 
voltage regulator and speed governor, vshich moans control 
of both terminal voltage and rotor speed and used the 
transfer function method of analysis, 

Q 

Root locus method^ was usod to study the effect of 
regulator parameters on the stability by plotting the loci 
of the roots of the system characteristic equation. A 
qualitative assesmont of the system transient performance 
can bo made by the study of the root locus diagrams. 



1 o 

Using the system sensitivity equations , the 
sensitivity of each of the charaoteristio roots to 
different parameter values can be obtained and a region 
in the parameter piano can be determined within which the 
system is stable. 

1 1 

Using generalized Hitrovic method formulas are 

obtained for the calculation of optimum loop gams and 

dan^ping gains of a certain clacs of voltage regulators 

for better steady state stability. The generalized 
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Mitrovic method is also employed to ootain a prescribed 
root configuration in the s-plane by vara at ion of the 
regulator parameters. 

In most of these methods tne system model is obtained 
as a single higher order diffomntial equation, including 
the regulator dynamics* Using cla.ssical control techniques 
the system characteristic equation is tested for stability. 
The parameters aia varied till a satisfactory performance 
IS obtained. In this chapter, the second method of Lyapunov 
v/hich relies heavily on the state variable formulation, is 
used to obtain optimum regulator parameters . The design 
criterion is that the settling time of the system variables 
should be minimum in the event of system disturbances. 

3*3 LYAPUNOV FUNCTION AND MINIMUM SJUTOLINO TIME 

The Lyapunov's second method can be used to design 
optimal regulators which ensure an ajsymptoticalDy stable 


system. Also the system output v/ill be continuously 
driven to^rards the desired value. The Lyapunov* s method 
requires the utilization of a continuous scalar function 
of state variables called the Lyapunov function in 

con 3 unction with the system stale equation. Depending 
upon the properties of V-function and its time derivative 
V(X), the stability or instability of the equilibrium 
state can be proved, 

A Lyapunov function gives a measure of the system 
slalo at any given instant of tme and hence it can be 
considorod as a measure of the distance of the state of 
the sysbom from the equilibrium state, in tho sbato space. 
If a V-fuiiction is loaown for an autonomous system, then it 
can bo used to estimate the rapidity of the ■cransient 
response or tho rate at which tho state comes bach to its 
normal state from tho disturbed state. A measure of the 
transient response can bo taken as tho normalized rato at 
whioh the V-function changes. 

Considering the origin as a stable equilibrium 
1 5 

state, let 



in some region of the state spaco excluding the origin. 
Equation (3.1) is intograbod between limits 0 to t assuming 
i\ to be constant in the specified region, to give 



(3-2) 


t 

V(Z) 7(2o) exp(- / ndt) 

0 

G?lio above equation gives a measure of how fast the origin 

is rc-ohod from any givon initial state The time 
14 

const, nt for •che variation of V-f unction can be considered 
as l/n. Thus a large vdlue of n corresponds to faster 
response. For a given system, a large number of V-functions 
can bo determined. Hence the largest value of n should be 
taken as a figure of merit for the system transient perfor- 
mance. 


For a linear time invariant system 

Z * X, X(o) * Xq (3.3) 

a Iiyapunov function can be detorimned easily as a positive 
dofiniLo quadratic form in the siito variables, with a 
posioivo definite Q matrix, as 

V(Z) « z'^’ P Z (3.4) 

where P IS the positive definite solution of 

E + A^E + « = 0 (3.5) 

and v;ith a negative definite Y-function as 

V(X) « - X® Q Z (3.6) 

The figure of merit for the linear system given by 
equation (3*3) can be defined as^"^ 

Tj «* Min(Z^ Q X) 

Z 


(3.7) 



Lo tho constraint 

= 1 (3.8) 

The mmimiaation of n can bo achieved by using Iiagrangian 
multii^lier technique, with the lagrangian multiplier X . 

The Hamiltonian function is formed as 

H(X, X)=r'^ QX+X(1 (3.9) 

Minunissntion of n is the same as minimizing equation (5,9) 
\;itVio iii -niy f'cnrstr -mt. Thus laiiaimizing H, gives 

(Q- XP) X = 0 at X s= (3.10) 

Houc'o 

1,0. n»X>0 (3 *12) 

Prom equation (3.10), x is shown^^ to be the eigenvalue 

^ -j 

of ><# P' . Thus ti can be taken as the minimum eigen value 
of ^ P**^ . Since a large number of Y-functions can be 
found for the system, n should be taken as the largest value 
from the set of minimum eigen vulues of Q P" . The problem 

is tho same as finding the minimuLi Vv-'luo of the largest 

-1 

eigen value j of P Q , Thus the figure of merit for the 

1 6 

linear syoLom is obtained as 


n = Min ( )] 


(3.13) 



3.4 SOIUCTOIT OP THE PROEDB5I 


Tho system shown in Pigure 2.1 is provided with 
simple voltage regulator and speed governor, !Dhe system 
model in linoar time invariant form as obtained m 
Chapter II is given by 

i = A X + B u (3.14) 

She volbage regulator and speed governor dynamics are 
choson as 

Voltage Regulator: 

»ra = (3.i5) 

Speed Governor: 

o 

The regulator equations (3.15) aid (3.16) are linearised 
about the same operating point considered in Chapter II and 
thed augmented with the system model given by equation(3.14) » 
taking as additional state variables* Thus 

the sLauo space model for the conbrolled system is obtained 
in tlio form of equation (3.3)« 

The problem facing solution is that the voltage 
rogulauor and speed govemor gams and Kg respectively, 
have to bo optimized so that the settling tine of system 
variables is minimum in tlfiO event of any system disturbance . 
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The timo constantoT^ and Tg are selected as T^ « 2,0 and 
s= 0,1 , Por the operating point considered in Chapter II 
with rot. or angle of 26.3°, the matrix is given in the 
next T)Qii,o * 

The Q-mntrix is chosen as -Q = dia(1 0,1 0,1 0,1 0,10, 

10,10,10,10), and the matrix P is calculated from the 

e qua cion (;,5). The optimum values of and Kg are 

obtained by minimizing the largest eigen v^ue of P <i • 

-1 

Bvoluabio.i of bho largest eigen value of P Q instead of 

findJiig tho smallest eigen value of Q P , avoids the 

-1 

comput j bion of P in each iteration and hence saves much 
compuuint: time, Tho optimum values obtained by Rosonbrock's 
hill climbing technique^ (Appendix A) aro given by 

ICy « 3.7584 Kg * 4.711 

and tho minimum veluo of tj » 170,19465. The system with 
tliooo opti'nal paramebors is simulated on the digital computer 
for tho branoiont response \/hon there is an distur- 

bance in Ino field current of 0,05 p.u. The minimum settling 
time response is shown in Figure 3,1* The superiority of 
this response is established by comparison with Figure 2.7* 
file sotbling time is very less compared to the previous one. 


Tho voltage regulator configuration with sbabilizer 
shown in i'agure 2.5 and governor of the linear type given 


by 

^i* Kg 


w - Wq 


(3*17) 






IS then considered for selection of optimal parameters. 

Ihe g-iiis Kg and are assumed to bo the variables of 
incc'Cfj'^. The other regulator p£a?ameters are taken as 
given in Oliapter II. The optimal gains are obtained as 
S5 w.O and » 273 *0. With those parameter values the 
minimum setbling time response is obtained as shovm in 
Figure 3 *2. Only simple regulator configurations as given 
by equal ions (3.15) to (3.17) arc considered to illustrate 
the priiiciple. This method is applicable even when more 
than two variables (adjustable gains and/or time constants) 
are involved in the minimization procedure. 

3.5 COIKJIiUSIOir 

The various methods of selecting conventional 
rogulafcor parameters are briefly discussed. The second 
method of Lyapunov is used to obtain the optimal regulator 
gains. The performance criterion is taken as the minimum 
settling time without much overshoot. An improved perfor- 
mance can be obtained by this design procedure. The 
stability of the system is ensured automatically by a 
positive definite solution of the lyapunov matrix equation. 
The method con be extended to nonlinear systems easily. 
But the dotermination of a suitable lyapunov function for 
such high order nonlinear systems is difficult. There 
are no systematic methods which can be applied easily to 
obtain 7-functions for such nonlinear systems. 
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In bins chapter the configuration for voltage 
rogulacors and speed governors are assumed apriori and 
indcpoiidontly and then their optimal parameters are 
obtained. Hither to bhis has been the practice. However 
an Integra bod form of control lav/ v/hich is a combination 
of all the state variables of idle system which provides 
an optimal response will be very useful because it does 
not require any aprion knowledge of the regulator 
configurations. If this control law is linear and time 
invariant then it will be possible for practical 
implement at lonc. The following chapters deal with the 
do berminabion of this type of controllers. 



CHAJIBR 17 

OPTElAIi COITTROL OP STICIIROIOUS MCHHaB 

4.1 IiniRODUCIIOH 

In this chapter a now approach to the design of 
esccibauion and primemover control of the synchronous machine 
IS considered. The optimal output regulator is discussed 
for both finite and infinite time interval of optimization. 
The optimal control law is obtained as a linear feedback 
control of the complete state vector, using the linear time 
invariant state model derived in Chapter II. The perfor- 
mance of the optimally controlled system is compared with 
the convonbional regulator response for impulse type 
disturbances. The dynamic performance of the system is 
also investigated for large disturbances such as line 
reclosuro. The optimal ooutrol law is obbaaned at different 
operating conditions and performances at different loads are 
then compared. 

4.2 POEtlDIiATOT OP 0PTBI4I COHTEOE, PROBLEM 

The feedback control of synchronous machines has 
boon given a strong impetus by the modern optimization 
theory as developed by Pontiryagin, Kalman etc. which 
relios heavily on the sbate space formulation. It is 
usually difficult to translate the given system specifi- 
er bions in the conventional design procedure. However, 
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in the modified approach, an optimal control law is 
obtained by a suitable choice of a performance criterion. 
Ono of the difficulties of controlling a nonlinear system 
with nonlinear controls is that the optimal control 
policies ore not generally easy to implement. 

Ihe system shown in Figure 2,1 is identified as 
an optimal control problem. The objective is to obtain 
a linear constant feedback coiitrol law which will trensfer 
the system from the given initial state to the desired 
state* In so doing the control pnrstem must satisfy the 
requirements relating to the performance of the sjrstem 
such as the decirod response , desired control oJ'f ort etc . 
and also its implenontation* Tho porfornmeo index is c 
mathematical model of the perforrimce requirements. Ilary 
performance indices have been proposed in the iitrr«-> . 
The quadratic porformanoo index in the outpub ixid ooiiurwl 
vectors minimizes the error m bho ouLput vpritWoo id 
the control effort. Also it results in a linonr conorol 
law which IS easy to implement on preciical systeris. 

21 

The performanoo criterion is thus taken as 
1 ^ f . T’ T ^ 

J a 4 / (y^ Q y + u-^ E u) dt (4.1) 

^ 0 

where Q a positive semidefinite matrix, is the weightage 
associated with the output variables and E a positive 
definite matrix, is the weight age associated with the 



concrol variables. A unique set of weighting factors to 
saciefy tho prescribed design specifications generally does 
noL o^ist* However, the lack of uniqueness of the v/eighting 
factoxo doos introduce a flexibility which makes selection 

f 

of the performance index simpler. Depending upon the 
relac:ve importance of the stace or output vector and the 
control variables, the Q and H matrices are chosen to 
rolloct rhe desired closed loop sjrstem performance. The 
cons or .'inis on the control and output variables can be 
indiroctly taken into account by assigning suitable 
ponalijuG to the constrained v.-'riables. The and R 
piaui’ioc.3 lijiomsolVGS can bo chosen to reflect these 
pens] bjoo. 

'I’ho problem posed here is tliat once the performance 
oritorion is selected, it is required to find an optimal 
control law whic h minimizes the porformance index subject 
to tho constraints 

Z a A Z + B u (4.2) 

I s: C Z (4.3) 

The optimal control u can bo obtained by the application 

of Pontryagin*s miniiaum principle or the Hamilton-Jacobi- 
■« 

Bellnan theory. Tho existence of the optimal control 
requires tho prior investigation of the system output 
controllability as discussed in Appendix B. The optimal 
control law is obtained ' (Appendix 0) as 
- R’’’^ B*^ P Z 


u *t 


(4.4) 


\flicr 0 P 3.8 tho solution of natrnjs B.icca'ci. eq,uatioii 

P + x* 1 + A® P - P B R"”* P + C® <i C = 0 (4.5) 

v.ith PCb^) = 0, Por iiifinito time rogulator problem, the 

solution of algebraic Rioca'ci ectuation 

P A + A® P - P B R*”' B^ P + <i C = 0 (4.6) 

gives bho control law, Ihe op'cimally continued system is 
bhen C3.von by 

X a (A ~ B R’**' B^ P)X (4,7) 

Even an Uxis table system can be stabilized with the appli- 

1 7 

cation of optimal control . 


4.3 SOLUTION OP RIOCiiTI EQUATIOI 

The solutions of Riccati differential and algebraic 
matrix oriuotions are discussed below. 


Infinite Iimc Problem: 

Por the solution of algebraic matrix Riccati 

1 Q 

equablon, the method of euccossivo approximation ^ is 
used. In this method, approximation in contreiL policy 
space IS combined with stability considerations from the 
second method of lyapunov, A sequence of suboptimal control 
functions are then generated which have a monotonic 
convergence. Ihe initial choice of P is made such that 
the controlled system given by equation (4*7) is stable. 

Ihc Riooati equation at the kth iteration is given by 


P^^ A^^ + (A^)® P^ + ^ 


(4.8) 


a 0 



where 


= A B R"”* (4. 9) 

and. 

^ 0^ Cl C + P^“’' B R"'* B® P^'^ (4.10) 

Starting, with an initial value for P, the iterations are 
carried out until the difference between P^”^ and P*^ 
elements is within say one percent . 

Pinito Time Regulator: 

The matrix Riccati differential equation is 
integrntod backwards in time from the final time t^, with 
P(t^) ss 0 till zero time. The mabricos arc stored at the 
different time insliOnts and used for the feedb'^ck control. 

It has boon shown'^^ that the solution matrix P approaches 
a steady state value for t^^ sufficiently large. This value 
of the steady state matrix P is used as a constant feedback 
of states for all times from t » 0 to t s= t^. Also the 
performance deviation by using this control is not very 
much from the exact control law^-^. If the time varying 
control law is used, then the implementation of the feedoack 
law becomes difficult and requires a preprogrammed control, 

4.4 EERFORMANOB WITH 0PT31I/\l! REGULATORS 

Por the system given in Chapter II, the performance 
index with the following wuightago maurices arc orosen. The 
wei^tage on the output variables a5 > AW and is selected 
as <il » dia(1 0,1 0,10) and the weightage matrix on txiS control 



variables *2^^ and is selecced as R = diaCl,!)* Ror 

the lixCijjate time regulator the Riccati matrix is solved 
by the method of successive ap^jroximation the B-matrix 
is obbainod as 


P s 


11,424 0.0280 0,6950-0,6570 0.6410-0.0560 0.0580 
0.028 0.0610 -0,0250 0.0280 -0,0250 0.0029 -0.0089 
0.695-0.0250 2,1400-1.9660 1.9720-0,0134 0.0194 

-0,657 0.0280 -1.9660 1,3100 -1.8100 0.0129 -0,0187 
0,641 -0.0250 1 .9720 -1 ,8100 1.8200 -0.0126 0.0185 
-0,056 0,0029 -0,0154 0.0129 -0.0126 0.0240 -0.0187 
0,058 -0.0089 0.0194 -0,0137 0,0185 -0,0187 0.0178 


1 


For the finite time regulator problem, the finite time is 
taken as tj^ « 5 seconds. The maurix differential equation 
IS solved baelcwards in time, giving a steady state solution 
as 


P 


11,456 0.0286-0.7340-0.6720 0.6760-0.0560 0,0590 

0.0286 0,0608 -0.0249 0.0283 -0.0251 0.0029 -0.0089 

0.734-0,0249 2.1414-1.9674 1.9728-0.0144 0.0208 
-0.672 0.0283 -1,9674 1 ,8122 -1 .8134 0.0139 -0.0199 
0.676 -0.0251 -1.9728 -1,8134 1.8179 -0.0135 0.0197 
-0,056 0.0029 -0.0144 0.0139 -0,0135 0.0240 -0,0187 
0,059 -0.0089 0,0208 -0.0199 0,0197 -0,0187 0.0178 


The controlled systom with iho optimal regulator is 
simulated on IBEI 7044 digital computer for the dynamic 
performance when these is an type disturbance in 



the field current. The fourth order Sunge-Kutta integration 
me'chod is used to solve the equations. The response of the 
optimally controlled system is shown m Figure 4.1 for the 
infinite time regulator and in Figure 4.2 for the finite 
time regulator. The co]i5>arison of Figures 4.1 and 4.2 with 
the conventional voltage regulator and speed governor 
response given in Figure 2.7 reveals the superiority of the 
optimal regulators. The response decays exponentially with 
no overshoot and oscillations and takes much smaller time 
to settle down. 

4 .5 EESPOHS® OF TI-IB SYSTM FOR LAR&B DISTTTHMHCBS 

The design of the optimal regulators for the single 
maoliine system is repeated at different operating conditions. 
The average value of the performance index is show/^^ to be 

J s tr (F) (4.11) 

The values of S are tabulated below for different rotor 
angles. 


*0 

17.700 

26.5® 

58.00 

77.00 

87.00 

— * 

j 

16.99 

17.3 

17.4 

16.77 

16.12 


From the above results, it can be concluded that the 
performance value with linear time invariant feedback 
control obtained for one particular operating condition 
almost remains optimal over different load conditions. 
This ensures the adaptability of this control even for 
large disturbances. To prove this point, the performance 






of the system with optimal regulators is investSgSPfiBTwhen 
the second transmission line in figure 2.1 is reclos^* The 
system is originally operating with one line in service 
having the following operating conditions : 

*0 = 28.7“ = 314.0 ^a<, = 2.56 = 0.91 

igo" Xjjjj, = 0 . 0.323 1 .465 

When the second line is reclosed the problem is to transfer 
the system state Z from X|(o) to ZgCt), The responses of 
the system for this state transfer with optimal regtilabois 
computed at the post disturbance steady state operating 
point and with conventional regulators discussed in 
Chapter II are shown in I'igures 4.5 and 4.4 respectively, 
ffron these responses, it is seen that the optimal regulator 
is better than conventional regulators. Thus the optimal 
regulator calculated at a particular operating point can be 
used for different load conditions rath less performance 
deterioration. 

4.6 COHOIiUSIOir 

Jin optimal output regulator for the system considered 
is obtained at different operating conditions. The 
superiority of the optimal integrated control is establi- 
shed by oomparison of responses of optimal and conventional 
regulators* The performance of the system is investagated 
for large disturbances* The optimal control law requires 
the entire state vector for the purpose of feedback. 





CHAKDER V 

A DYITMIC OBSERVER EQR STNCHRONOUS WAfirTTi^iTi) 

5,1 OTROBUOTIOR 

The optimal regulator designed in the last chapter 
calls for the direct measurement of entire state vector. 
This is impractical in many problems because the state 
variables chosen need not necessarily correspond to 
physic.'lly measurable quantities. In some oases the 
motiourement of some of the variables may be forbidden. 

Even if it is possible^ it may not be economical to do so 
in many cases • As mentioned earlier it may not be possible 
to know bhoir post fault steady state values to compute the 
deviations* Therefore it is necessary to reconstruct the 
staues by employing either a Kalman-Bucy filter or a 
lueiibercer type observer, from tlae meesurements of only a 
few output variables. 

In this chapter a dynamlo observer of the Luenberger 
type is discussed which reconstructs the state vector from 
the available output measurements. The performance of the 
optimally controlled system obtained in Chapter IV, cascaded 
with the compatible observer is then obtained. The transfer 
function matrix relating the input and output of the system 
is also derived. 
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It is shown that for a linear, time invariant, 
finite dimensional I dynamic observable system, it is 
always possible to obtain a compatible observer* A 
compatible observer is one whose output equals the 
stale of b]ie system to within an exponentially decaying 
error * 

5*2 DmMIIC OBSERVER 

The problem is that given the linear time 
invariant oystem 

3C as A X 4* B u (5*1) 

Y » C X (5.2) 

it is required to recons uruot the unavailable state 

variable 0 , The system has to be both output controllable 

and observable as discussed in Appendix B* Eor the 

linear system described by equations (5*1) and (5*2), the 

22 

Observer ils defined as 

Z » PZ + GY + Hu (5.5) 

1 » Y + li2 Z (5.4) 

where X is the reconstructed state vector. The observer 

state Z is p-dimensional where p » n>^m, n being the 

system order and m being the number of available outputs. 
22 

Por a choice of 

H » T B (5.5) 

where T is the solution of 

s G 0 


T A - P T 


(5.6) 



iho otsorvor and syateu states are related Tjy 




Z a 2 X + e 


(5.7) 


where o lo the error vocbor bet\/eon the observed and actual 
stale vcobore and is given by 

e » Qxp (r b) e(o) (5.8) 


and e(o) lo the error at at t * o. The reconstructed 
state IS given by 

X a L(liiQC+ ©) (5.9) 

where e is defined as 



with Ojjj ao an ai-dimensional null vector, 
partitioned natrices ei’ven by 


( 5 . 10 ) 

1 and W are the 


L a (L^ I Xg) 


(5.11) 


W 



( 5 . 12 ) 


In the J^ovB treatment it is assumed that S' is a 

stable matrix and S' does not have any eigen value in common 

with that of system matrix A. In order to completely 

specidfy the observer dynamics all the above discussed 

matrices arc to be determined. The matrices I and S are 

chosen such that the pair (I?, G) is controllable. The 

22 

matrte T is obtained in a straight forward manner as 



58 


It = -r''(r)nns (5.13) 

where 

R « (G ; r G : G : g) (5.U) 

xs tlio oorrbroHabxlxty matriz of the observer system^f 


= Lo’ 




... (5.15) 


xs the observability matrxx of the system given by 
equations (5.1) and (5.2), 


** i 

^ ^ “ xio ^ “n ® ^ (5.16) 

xs tho oliaraotoristxc polyromial of the system matrxx A 
and etj^’s nro tho coofficxento of its characteristic 
equation and 


n 


*1 ^ ”2 ^ ®n-1 \ 

*2 ^ “3 ^ “n ^ 

®ft-1 ^ *n^ °m 

*n ^ ®m ••••* ®m 



and Ojjj are mxm identity and null matrices respectively, 
22 

It has been shown that the observer is compatible if the 


matrix 7 is nonsingular. 



5,5 S5niCir.l01T0US machihe observer 


For the system shovm m ligure 2,1, it is required 
to reconstruct the seven state variables from three output 
variables namely the rotor angle, rotor speed and machine 
terminal voltage. As discussed in the last section, the 
observer order becomes four, Eor the operating point 
considci-'ed in Chapter II with rotor angle of 26.30, the 
matrices A, B and C are already given. The system is 
found to 00 both output controllable and observable using 
equations (J#4) and (B,5). 


Tho observer matrices E and G are selected as^^ 


F « 



0 0 0 

-20 0 0 

0 -20 0 

0 0 -20 


Qt rr 


1 

0 

0 

1 


0 

1 

0 

0 



The matrices jZf(P), n and S are calculated using 

equations (5*14) to (5*17)* Then the matrix T is obtained 

from equation (5*13) and the W-raatrix is constructed using 

tho matrioos 0 and T, The submatrices l-j and I12 are 

22 

obtained by assuming that B ss W * 


The H-matrix is 



oalculal/Gd using ociuation (5«5)» Thus all the observer 
mstncos aro obtained and aro given by 




0*0484 -0.0025 
0.0315 0.0491 
0.0079 -0.0027 
0,0079 -0.0027 


-0.0061 0.0052 
0,1226 -0.1048 
0.0268 -0.0249 
0.0268 -0.0249 


-0.0049 -0.0119 
0.0989 0.2376 
-0.0234 -0.0068 
-0.0234 -0.0068 


0.0159 

-0.3189 

0.0116 

0.0116 


H 


“*-0.0029 

-0.1277 


0.0597 

2.5550 


0.0093 

-0.0139 


0.0093 

-0.0139 


T 



0 4.2175 526.90 264.90 4231 .60 3063.55 

1 -0.0068 -0.0400 -0.0689 -0.3218 -0.0968 

0 0.3068 4,7506 3.8470 11.558 8.3789 


« 


0 0 
0 0 
0 0 
0 0 




-86.52 

-10530.0 -5289,0 

-3 .496B9 

-61270 

-4.503 

-525 .40 

-262.40 

-84630 

-3064 

-2.896E9 

2 .801E8 

3.496B9 

-4224 

1 .54B9 

2 .896E9 

-2.801B8 -3.496E9 

4224 

-1 .54E9 


In tho niC''lrlx E'*’ ifloans 10*. Once f^ese matrices are 
detonnincd the observer dynamics aio completely specified. 


5.4 PERrOTITlUTCE Him DmillO OBSERVER 

The observer constructed in the previous section is 
cascaded wluh the optimal infinite time regulator determined 
in Chapter IV. Tho response of the system for a disturbance 
of 0.05 P.u. in the field current is obtained by simulating 



the ca&cadod system. The response is shown in Figure 5,1, 
The responses have more overshoot in the initial portions 
but decays fast to the steady state values. 

The observer matrices are found^^ to be sensitive 
to the system operating conditions. Hence for each 
operating point f the observer dynamics are different and 
therefore the control law is different for the various 
operating ooiiditions , But this is not the case with all 
the olf'-uo vni'iablos available for iijeasurements as 
discussed in Luo last chapter. The error in the 
r 0 con siU'U clod state vector decays exponentially depend- 
ing upon Lho values of the elements of F-matrix as given 
by equauion (5*10), By choosing large negative eigen 
values for F, the error can be reduced; but this 
introduces largo goins in the transfer matrix between 
input and output. By increasing the v/eightages on the 
output variables the system can be brought to equilibrium 
state in n quicJccr time at tho expense of large inputs. 
Henco a oompromiso must be sought between these two, 

5,5 TR/'drSPlE T-I/lTBIX 

TUo braasfor natruc^^ relating the output and 
input arc derived by assuming that all the initial 
conditions arc zero. The various equations are written 
in Iiaplaoo transformed variables as follows : 

From equation (5*3) 

Z(b) SB (si « F)"** [ G Y(s) + H ^(s)3 


(5.18) 




From equation (5*4) 


X(s)»Ii'jY(s) + 3 j 2 Z(s) (5*19) 

From oquatioii (4 •4-) 

U(s) = r XCs) (5.20) 

where 

F « - R P (5.21) 

Assuming tliac X(s) * X(a), the following equation is 
obtained 

U(s) « F[Li YCs) + I 2 Z(s)] (5.22) 


Substituti*;, foi** Z(s) from equation (5.18), the equation 
(5.22) rcUuco*? to 

U(8) a ' 111 + I 2 K(s) G] Y(s) + F I 2 K(s) H U(s) 

... (5.23) 

where 

K(8) - (si - P)"’ (5.24) 

From equation (5.23)» U(s) is obtained as 

U(s) » Cl - F I 2 K(s) hH + ig K(s) Gjy(s) 

... (5.25) 

IhuB the transfer function matrix ]S'(s) in the relation 
UCs) » H(s) y(s) (5.26) 


for the operating conditions considered, 
equation (5.25) as 


lJ(s)gi" ^— " 

DCs) 


T45.5e^+2205s+38 0.02s^-s-2096 
355 8^+71400-777 ( -3 .2s^-462 .5s 

-7311) 


IS obtained using 

0.27s^+9.5s+16.5"’ 

0.63s^+73b+1220 



whero 

D(s) « + 57.85s + 354.32 

Thus the i,jth olement in the transfer matrix N(s) refers to 
the tratisftsr function bot’veon the jth output and ith input 
and this enn be practically implemontod , However, this 
control ^ives optimal response only for the specified 
operating point. The transfer function matrix H(s) is 
d afferent for the different operating conditions. 

5,6 conclusion 

oonpatible dynamic observer is obtained for the 
synchronous machine system bo reconstrucb the unavailable 
slatou* Tho cffocl of observer dynamics is that an 
exponentially dto.^ying error is introduced in the 
estimated value of the states. The observer dynamics 
arc very muoh sensitive bo the operating conditions .Thus 
for onoh operating point the observer dynamics have to 
be modified. The input output transfer functions are 
dorivod \/hich indicates tnat time varying functions are 
introduced in the feedback paths. The gams of these 
feedback pabhs are different for different operating 
conditions. The response of the cascaded system is 
obtained and compared with complete state feedback 
system rosponso* 

Thus oven though the difficulty in the measure- 
ment of post fault steady state vector is avoided by a 



choice or measurable output variables whose steady state 
values Tor all operating conditions can be obtained apriori 
the majoiT defoob oi the control scheme proposed in this 
ohaptoT^ LG b’ut tho transfer function matrix is difforenb 
for difiorent operating conditions. This makes the optimal 
control law useless for large perturbations unless a pre- 
programnod controller is used. This of course will be 
practically difficult to implement. To overcome this 
difficully suboplimal controla which give rise to 
rcaoonailj'' i'ood. rcoponcG and also \;hich are linear and 
tine Lii/. i’Lfmb, arc considered in the following chapter. 



CHAPTER VI 

SIDOPTH'IAL OOKOEOl OB’ SIHCHROirOirS MAOHUSfE 
6.1 IHTRODUOTIOT'; 

In rooont years j the conecpt of state space and the 
use of Po-ioryasin’s minimum principle have resulted in 
analytical techniques for the synthesis of optimal regula- 
tors for nullivariabOB systems. The resulting optimal 
control law c.ills for a complete measurement of the state 
vector. The state variables being mathematical variables 
introduced Cor the convenience need not necessarily 
correspond to physically measurable quantities. Thus one 
must consjdor the problem of controlling a physical system 
with 03ily some of Lho state variables available for 
feedback. 

Often it is sugjscsitod that the unavailable state 
variables Ovin be rcoonatruotod via a Kalman-3ucy filter 
or a state Roconstructer^^. But this introduces transfer 
functions^® in tho feedback paths. It is quite impractical 
to reconstruct tho state variables m large multivariaole 
systems such as those which occur in interconnected power 
systems f ohomloal process control etc. Purther it is 
shown that the resulting gains and time constants in the 
feedback path depend on the operating conditions • Thus 
it would bo better if the control variables are chosen by 
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8 , linear combina'tioii of "bhe available oubpiit variables 
instead of employing an observer^ This may result in 
detenora oion of the porformance of the system* However 
if Iho perfOiiui^oe is reasonably close to the optimal 
one then lu \o.ll be better to go in for this linear time 
invariant; conlrol law* 

In lihis clidptor, a suboptimal control technique 
is outlined \/*ioro in it is not necessary to reconstruct 
the unnvt' liable state vaiiables* Thus the merit of this 
ap.'troaoh lies in the fact that a linear time invariant 
control la\/ is obtained ivhich is easy to implement. The 
performance oC the system with suboptimal control is 
dolorminod. The feasibility of using a particular control 
law over \ado range of operating conditions is then 
discussed • 

6,2 STATDIiU^L' (T THE PHOBIiBM 

The problem posed here is the following: G-iven a 
linear tiiiio invariant dynamic system 

X A X + B u (6.1) 

Y a.CX (^*2) 

it is required to find an optimal control law which 
minimizes the quadratic performance index* 

J*i /**(2?QX + u^R u)dt 
^ 0 


(■6.3) 



In addition, the control law is constrained to be a linear 
function of oho output variables as 


u 


3? Y 


(6.4) 


If all tho Ooito variables arc a'v'ailablo for feedback 

n-l 

( 1 , 0 . if o’" exists) then the optimal control law for the 

1 7 

state rogulabor discussed above is given by 

u « r X 
where 


1* 




J I 




(6.5) 

( 6 . 6 ) 


and P JS tlio pooitivo definite solution of 
PA + '^P-PBTl"‘'’3®P + a«sO 


(6.7) 


JT 


The optimal v iluo of bhe porformance index is shown to 
be 


J s X(o)’^ P X(o) 


( 6 , 8 ) 


Subopt imcO. Control* 

Phe control luw is constrained to be a function of 
output variables, Then using equation (6.2), the control 
law becomes 

u « P 0 X 


Then the porformaiice index is given by 

J ■ y /* r'(9 + / E t o)z at 

and the closed loop system becomes 
X a (; + B P 0)X, X(o) » Xq 


(6 JO) 


( 6 . 11 ) 
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UoWj Lho pPobl©n is “to ilmd tlie Glemsnlis of "blie ccaa.'trol 
matrix P, \/hioh minimizes the performance index J. Sub- 
stitutin,', bho solution of state vector obtained from 
equaoion (0.11 ) in equation (6*10), the performance 
is obb. mod v'ls 

J s. ^ X(o)^{ / jZf^(t) [Q + 0® R p C ] jZf(t)dt }X(o) 

( 6 . 12 ) 

where 

0{t) a oxp (A + B P C)t (6.13) 

is the ruiiaat.iont£il intrix for tiio system 2 >^reA by 

equation (6.11). Prom equation (6.12), it is seen that 

tho porronianoc Index is a function of both the control 

matrix P and tho initial state Z(o), To avoid the 

dopendenoo of J on X(o), the initial state o'-in be ireated 

as a random vector uniformly distributed over the surface 
2b 27 

of unit Gpiioro * ‘ . Then tho 'vor^ge vetLue of the 
porlormance ii^ox is given by^^ 

^ “ 2n ^ (Q + 0® P® R P 0) j2^(t)3 dt (6.14) 

Thus the problem reduces to tho determination of a sub- 
optimal control matrix P, which minimizes J subooct to tiB 
00 ns tra int e quat ion ( 6 . 1 1 ) . 

Using the trace minimization procedure, the solution 

26 

of tho subop bimal control problem is shown to be 

s = - n"’ b''' K 1 0® (0 1 0®)"'' 


(6.15) 
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where K is vuo positive semidoxinite solution of 

K + *'1 K + Q + C® E P C = 0 (6.16) 

I is Uio pooitivo dofinibo solution of 

II a| l- Ii + I as 0 (6.17) 

and A^ is a stable natrix ^ivon by 

=* (*- + B P C) (6.18) 

iUlio subopt iaaL matrix P is obt^ mod by an iterative 
al(?orithm'* usiiig oquabions (6,15) to (6,18). An initial 
Vviluo of P is ehooon such tht'it the matrix A^ is stable. 
Then naLrious K a d L arc solved and bhen a now value of 
P 10 obiuinod using equation (6.15). 3Jho process is 
repeated until bho difforenco between successive values of 
the olcnonbo of the rartrix P is within one percent, 

ir 0 oxisbs, then bho control law obtained is 

the same cti ICalraan’s optimal regulator, The expected 

26 

value or tho performance mdex can be shown to be 

i » tr(K), (6.19) 

For the cowploto state feedback oasG» bho average value 
of porforrvijiioe index is 

J a tr(P) (6.20) 

26 


It Is also ohown that' 
tr(P) < tr(K) 


( 6 . 21 ) 



Heno© tho conbrol law oblniaGd using ociuation (6*8) is noi; 
optuiol but lb .all bo used as a suboptimal control law. 


6.3 SUBOPTU-tH. i23Gmi.IT OH PC® THE SYSTM 

Tho xoliowing operating conditions are obtained when 
the maohine is delivering rated KVA at unity power factor 
to the lAiTiiiito bus, using the phasor diagrEun given in 
Pigure 2,3 s 


6o » 

16,0^ 

O 

II 

• 

o 

^fdo * 

1.82 1^^ = 0.72 

^qo 

a 0.694 

^0 “ ° 

%cqo ” 

0 

V 

mo 

* 1 .092 

» 0.5« 

V 

qo 

a 0,945 


Tho syaton ’u. trices E and C arc calculated for the 

above oonditioao and aro given by 


0.0 

1.000 

0,00 

0.00 

0,00 

0.00 

o.oo" 


0.0 

-0,167 

-56,14 

14.56 

-36.14 

-79.56 

22.50 



1.067 

-2.13 

-36.23 

-24.15 

664 .23 

-362.30' 


-879.9 

2.240 

-0,66 

-76.08 

-12*07 

1395.00 

-760,85 


-419.0 

1.067 

1,33 

-36.23 

-3C.64 

664,23 

-362.30 


696.0 

2.720 

897.14 

-1345.70 

897.14 

-53,83 

-26.91 


^484.5 

2.040 

672.86 

^009. 28 

672.86 

-40.37 

-35 *88 


fo.o 

0,0 

2.' 

125 0 

,664 

-1 .33 

0,0 0, 

s 

1 0.0 
u. 

52.08 

0,< 

00 0 

.00 

0.00 

0.0 0. 

0 
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1 .00 

0.00 

0.00 

0.00 

0,00 

0.00 

0.00“' 

c » 

0.00 

0 

0 

0.00 

0.00 

0.00 

0.00 

0.00 


J-C.275 

0,00 

0,00 

0.285 

0,00 

-0.107 

0 , 00 _^ 


The liiatricus Q and R Xor iih© scate o^Tid coirbrol 

variables roopuctivoly are chosen as identity matrices and 

ulio optimal foodback fiiccati matrix for the infinite time 

state rcf.ulaoor is obt lined by the method of successive 

28 

approximaliun as 


? « 


2,270 0.009 1,5800 -1 .4500 1,4500 - 0.078 0.084 

0,009 0.079 -0,0016 0.0035 -0,0035 -0,003 -0,006 

1 , 980 - 0.0016 4 . 3800 - 4,0000 3 . 9900 - 0.009 0.015 

- 1 .^ 5 *' 0.0035-4,0000 3.6800 - 3.6600 0.071 - 0.016 

1 , 4 j 0 - 0.0035 3.9900 - 3.6600 3,6700 - 0.011 0.016 

- 0,078 0.0030 -0.0090 0.0110 - 0.0110 0.082 - 0.078 

1 0 , 0 CJ 4 - 0.0060 0.0150 - 0.0160 0.0160 - 0.078 0.087 


The subopt iioal control matrix P is obtained by the iterative 
procedure using equations ( 6 , 15 ) to ( 6 , 18 ) for the above 
opor^tinfs, conditions as 


F » 


- 0.5857 


— 0.5000 


1 .4761 



0,9353 


1.1428 


0.2071 



and 

the TC r, 

latrix is given by 





’"m .240 

0,014 

1 .420 

-1 .310 

1 .300 - 

•0.078 

0.081 


0.014 

0.022 

-0.002 

0.004 - 

•0.004 

0.002 

-0.005 


1 

-('.002 

4 .48 

—4 .090 

CO 

0 

• 

0.028 

-0.009 

K s: 

-I 

0.004 

1 

• 

o 

o 

3.760 . 

-3.740 

0.005 

0.057 


i .3('0 

-0.0i')4 

4.080 

-3.740 

3.750 

0.005 

-0.006 


-0.0711 

.'.002 

0.008 

0.005 

0.005 

0.0B7 

-0.081 


O.U0.1 

•rfWI 

-O.005 

-0,009 

0.057 • 

-0.006 - 

-0.081 

0.090 

My »' 

mnjJ'Lni. 

i iLa-' .1 

pfWiiLivo 

defiiiitc 

solutini 

’ for tho 

Ii-n.. 

iL.i’i.: 

loi.u; »; 

;c}i sL..tgu 

oi conputation, 

the 35’‘ot(jm 


ijL-ili.il i i..v L,i .•{Jij'Ucucl vvibJi Lli', juboptirial control. 


i >.4 pjijj.u'' 0 . 4 . u-vJaii wii’ji auBoraniffi oai'KOEL 

‘i’lio ,»v;4M'orjiianoo oi‘ the rjiuchino oystoH v/ith 

till. jjuhc-i'LUuU. J\jodback cart-rol law ia obtained for u]3e 
opor-' L hi.'. palnL coiioidercd. The rosijonoo o’;; the sub-- 
C'pL p.ciit!:;‘ol fjyol-CTn for -an diaturbanco i.n the 

riti-UI cui'-v.'i: of 0,05 P»u.. is aliovm in Pisurc fi-1 and tho 
reai'oi}:!< . Ith cojuplotc at.itc fccdbacl: for tlio c'rc 
ooml i L i.i-jjjj ia »!u)VM in I'’iguro G.2, Xho 'ivcivifo pcrforriaijce 
index V* liuj ;i.« j.-ioru only by 0.23 over t}:-: 'liyi-ri.:-! valuo, 

Tho i>e- ’i*< '’ti- notj deviation, by the uce of suooptr...’.iil control 
.i.rj ve ry lojin. Tlic cuboptirnl dcoign’is r^,pcatod for 
diffpront opor..tine conditions. The suboptimal control 
1 .V/ obtc.ined f«r one operating point is used tit different’ 
Oj[)i)rab±rvi oonditious and the average porforraunce index 
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FIG- 6*2 ReSPONSE WITH pPjiUAL 
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values arc calculated and arc given below 


1 

1 

i 

Expected value of J 1 

Optinril 

Suboptimal 

Suboptimal 
with P at 46° 

17.7° 

U.73 

14.91 

15.32 

26.3° 

14.52 

14.54 

14.64 

46.0° 

14.20 

14.45 

14.43 

1 


Fror.i the rcstiLts it can be concluded that control 

riatrix P c. IcuJ.atcd at ^ » 46° can bo used over a range 
oi' load ooiii. .1. Ijions between say 6 = 15° to 6 = 60° without 
much iiuLiVuTk -ICO dutorioration. This enoures that the 
subopbiu.l coabrol lav/ given by equation (6.9) is 
ro*iyonal)lj' ;;ood oven i’or large disturbances. 

■j’hij fj.ytil.oju povi'orHnnco for a largo disturbance 
duo iivj lino roclosure for tLc syatom shown in 
Pigiiro in disousuod with tho suboptinal control law 
obtained at Ihc post cUsturbanco steady state conditions. 
I’hc pvod'Jotivbancc opoi'ating conditions ai*e : 

=. 51 .(!« Wo • %o • 2-°°^ ^do " 

iqo"- "-52 ikao- 0 ^50 - ° ’'do “ ’qp = ■■ 

The jicrPorrianco curves ore plotted in Pigurc 6.3 for 
this aysten state transfer with suboptimal regulator 
and iu jL-'iguro 6.4 for tho sane conditions with conplete 


state feedback. 








In tac above analysis the output variables used 
for the feedbaok are A6 , Aw and ^ Y^, The steady state 
values of Aw and l^Y^ arc zero irrespective of the 
operating conditions and hence there is no difficulty in 
obtaining the doviw-\tions required for feedback. It is 
necessary to know the post fault steady state power angle 
to obt'-in the deviations of fi . This will be difficult to 
know opriori in large systems . However , if the machine 
operates nb constant load angle and the internal voltage B^is 
varued to deliver dirfcroat powers (as in the case of d.v 7 .r. 
uaohinoo which will bo dto cussed in the next chapter) this 
dil'iculty can be oircmx/enood . As an alternative the output 
varinbloo can bo chosen to bo A v/ and and the sub- 

optiml control can be obtai:ied. In this case the measure- 
ment of deviations of all the output variables will oe easy 
because their steady state values are zero for all operating 
conditions • 


I*or tho operating point discussed earlier when the 
nachine is delivering rated KVA at unity power factor to 
tho infinite bus, with both lines in service, the suboptimal 
control law in obuainod in ternB of the now sot of output 
Vvird ablos • Tho f feedback control matrix for this case is 
given by 


F 


5.658:3 -9.3192 

-27.3824 25.714^ 
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I’lie volue of tno perforuenoe index using equation (6.20) 

IS obuainod as 19*3941 foi* this suboptunal control and 
Llio deviation from the opLiual one is 4.8714. 

6.3 Ca^OLUSION 

ITho suboptimal control of the synchronous machine 
system IS discussed. The control law is constrained to 
bo a linear combination of the output variables, vhich 
arc avniLnblo for dirccl measurements. The feasibility 
Oi*' UG3i:g a particular suboptiml control law at different 
(jj>o I wing points is invos bigoted . The performance of the 
system with tne suboptii'itil control is compared with that 
of a cOiipQwtc feedback case. The performance of nhe 
ojstora \/ith suboptiiinl contx'ol is studied for large system 
disc urbane os. Tho implcmcntabion of the suboptimal control 
la\/ is siuplo as it requires tho measuremem; of only the 
avoilablo output Vviriablos. 



Cr\I>a)ER VII 


SfAlE SPACE MODEL OP 

DIVIDED wiudhtg- hotoh synchronous machine 

* 

7.1 IN'i'RODUOIION 

In this chapter, a d.w.r. synchronous machine with 
two field windings displaced by 60oE is considered, The 
staoe space model is derived in a form most suitable for 
the application of optimal regulator design. The response 
of the uncontrolled system is compared with the system 
response provided with angle and voltage regulators on the 
field windings and a conventional speed governor for the 
control of input torq.ue. 

With the growth of large na'cional "rids, power is 
transmitted in bulk ovox' long disbances at very high 
voltages. In rocenl times the use ox high voltage caoles 
is also in the increase. Thus the system has to supply 
large reactive power. Por economic reasons, modern 
generators are of lax go capacity and arc designed with 
low short circuit ratios. All these factors add to the 
problem of maintaining stability ospeoaally under leading 
power factor operations. A number of methods were 
suggested to overcome tne stability problem. One method 
suggested is tho use of synchronous machines with an 
additional field winding and suitable excitation control 
systems. 
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Sopper aiad Pagg*'^ considored a machine with two 

field windings displaced by 60®E, having ang i r regulator 

on one winding and voltage regulator on the other winding* 
31 

Kapoor ot*al* suialysed a machine with one winding on the 

direct axis having fixed control tjnt? , another winding on 

the quadrature axis with angle regulator, lamamurthy 
32 

et .al . used a machine \/ith two field windings, one in 

d— axis and another on the q— axis, provided with angle 

voltage rogula uors respectively on them. Krause and 
33 

'fowls conuidored a machine with two field windings for 
the synchronous machine danping. 

By all these studies, it has been establishod 
that such a machine has greater steady state and transient 
stability limits than the normal machine employing conven- 
tional regulators* fho limit of roacti've power capacity 

31 

for normal machines at no load can be shovm to be 
^m/*q* d-axis regulation has the effect of reducing 

x^ -oly. But the q-axis roculation modifies x^^ and thereby 
extends the reactive limit* Even under loaded conditions 
tho ar'aotivo limit is more for bheso machines. 

In the conventional synchronous machine, tho 
excitation phasor is rigidly fixed to the rotor structure 
and displaced by the roior angle from the synchronously 
revolving voltage phasor of bhe power system during steady 
state operating conditions. By having a second field 
\/inding the excitation pliasor can be freed from the rotor 
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structure and this pcjrmitG the rotor to occupy any other 
suitable posibion depending upon the relative magnitudes 
of the two field v/inding voltages. 

7.2 STATE SP*''0E MODEL 

The schematic diagram of d.v/.r. synchronous machine 
IS shown in Figure 7.1 and the equivalent d,q machine in 
Figure 7.2. The d.w.r. machine is connected to the infinite 
bus through a double circuit transmission line as shown in 
Fijuro 7 .3 . Tho general nonlinear system equations describ- 
ing the performance of the system aro given by : 


Direct Axis Flux Mnkageo: 


* *atd H *ard ^r ^ad ^d *" ^d ^d 
* kd “ *a1,d H + *ara Sr “ =^ad ^d * *liia ^ 


(7.1) 

(7.2) 


Quadrature Axis Flux Linkages: 

“ *atq H *" *arq ^r * *q ^q ■** ^aq ^cq 
^kq * *atq H " *arq ^ *“ ^aq ^q ^kkq ^kq * 

Torque Winding Flux Linkages: 

’*'t * *t H ^r ^r ^ ^atd ^d “ ^atq ^q ^atd ^kd 

*atq ^kq 


(7.3) 

(7.4) 


(7.5) 


Reactive Winding Flux Linkages. 


'r * x-tr ^ ^r “ *ard ^d ^arq ^q ^ard ^kd 


‘•arq ^kq 


(7.6) 
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Direct Axis Voltages: 

^ P ♦d - 'a ^ ♦q 

P *lca '*' ^Ted 

Quadrature Axis Voltages; 

^ ^ P \ a 

® ■ 5?^ P Plea ■'' pfea ^iDi 


I'orquo \/ Hiding Voltugo: 

V1, = P n It 

lioactivc Winding Voltage : 


‘'^’r “ Wq ^ 'I'r **■ ^ 


D qua lion of notion of Rotor: 


v/horc 


T- = i_ - ’I'- 2 .. 
e a q q a 


is the oloctrical torque in the air gap. 


troncratox' Torminal Voltage Condition; 



+ V 


2 

<1 


Transmission Line Equations: 

“ ^0 ""e h - ^0 ^q 

Vq_ « Cos 6+ rg i^ + i^ 


(7.7) 

(7.8) 

(7.9) 

( 7 . 10 ) 

(7.11) 

( 7 . 12 ) 

(7.13) 

^(7.14) 

(7.15) 

(7.16) 

(7.17) 
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Substxtutii^ ior tho flux 

(7.1) to (7.6) la the other perfonnanoe equations, they 
can bo put in tho vector Torm as 





E X 

a f(Z) + D U 










(7.18) 

wlioro tho state vector X 

IS chosen as 



X = 

(«» W, 1^, 


^q* ^kq 


(7.19) 

and 

the 

control vector u 

IS Civen by 



U s 

(vt» Vy., T 

sO? 

i' 



(7.20) 

TIto matrioos E and 

E nrn 

given by 




'l 

0 0 

0 

0 0 

0 

0 



0 

1 0 

0 

0 0 

0 

0 



0 

0 

*tr 

“*atd "^atq 

^atd 

^c*tq 


E a 

0 

0 *tr 


"^ard *arq 

^ard 

"^arq 1 


i 

0 

° *atd 

*ard 

-Xd 0 

^ad 

0 



0 

° *«tq -*arq 

0 -x 

Q. 

0 

2:a„ 

aq 



0 

° *ata 

*ard ‘ 

-Sad ° 

^IsOcd 

0 



0 

° *atq - 

*arq 

° -*aq 

0 

^q 








-i 







T 



0 

0 vfg 

0 

0 0 0 

0l 


E a 

0 

0 0 

% 

0 0 0 

0 j 



0 

m 

H 0 

0 

0 0 0 

oj 



Th" column vector f(X) is given as 
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w - 

-*o H 
-Wo ''r S? 

Wo ( Wa >+"(*81; 

Wo(Wg+r^iq) + wCxjiij-Xg^iji^-ii^i^-i^jij^j) 

-*0 ^ ^ca 

j;®o ^k, ijci 


Linear State Hodol; 

■ '^o diocuoood earlier when deriving a state space 
model for the normal syi ohronous machine, it will be 
convenient to linoarizc the uy3i<©m equations co implement 
the control Ihcorotic coxio-pts for -chc do&igii of opoimol 
rogulatora. It v/ill he jl'own later tnat the control law 
obtained will bo applicable oven for largo disturbances, 
itQ« when thu system is described by its nonlinear 
equations, TiioroXoru linearization of the system 
oqusition (7*18) about an operating point, loads to 

BXaPX + Lu (7.21) 

whore X and u oro the deviations of the steto and control 
vectors about the operating point rospectivcly. The 
operating point matrix P is given by 
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who 1*0 

^23 = (*atq^do’*ata^qo5/^‘ 

*24 = -^*ara*qo«i»rq*ao>'^ 

£35 “ t(X|j”*q)iqo*’*atq^o'*’^l*lC 9 o”*cirq*ro^ ^ 

*26 = C(“a-*q5*ao-=^*iEao-*Qta^to-''ara^ro]/“ 
Tgl - WqV^ Oos «£. 

*52 * *atq*to“*urq*ro”^q*qo'''*aq^lcqo 
*56 ■ ■"o^*o'“q^ 

*61 “ “Vo ®“ *0 

£52 = *a*ao"*ata*to"*ara^ro"*aa^o 


*65 “ "o^V^a^ 



fromiabiplio .tion of cqumon (7.21) by tho iivoroo of E, 
^TVDS the ^-snLom otaU rqu tion is 
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X a. A X H n u 


( 7 . 22 ) 


Ttio output V'Pi'.bloa ire oaoseii ?s bl:o rotor aji^e, rotor 
epood and tac machine terminal voltage . By linearizing 
oqurLnon (7.15) ^ivos 



'loo linourioing the oqun lions (7.16) e r\d (7 .17)7 

iVj » Vj, boo 4i!+ Aij - Xg 4i^ (7.24) 

^'’^4 ^ ‘Su* «o ■'^e + *0 (7.25) 


Thui’ ioro ooinbinlnc l-l*- ove tSirou cqut urns, 

4v„ - 0, 4« + 05 iij + 05 (7.26) 

^dlo^o 

“I • 7o(vao Oos Sui S^)A„o 

®5 “ (’^do + \o 

®6 - (V >^9 - '^do *o)Aa<, 


Thus tho output cquition for tho system can be obuainod as 


Y = C X 
whore 


(7.27) 


C e 



0 

1 

0 


0 

0 

0 


0 0 0 0 

0 0 0 0 

0 O5 Cg 0 


0 

0 

0 


f 

J 
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and 

y a (A5 , A w, 

Onc@ 1j1io operating poin't is soloctody the system matrices 
and C can bo ImmodDatoly cnlculabed. G?he complete 
system obate space model is obbamed as a linear time 
invariant model and is given by equations ( 7 . 22 ) and ( 7 . 27 ). 

7,3 PERPORIUIICE OP (CUE mTCOITTROILSD SYSTEM 

The system shown m Piguro 7.3 lias tiie follov/ing 
pardmotcrsi ''^I’oroin all uho quantities are in p.u. except 
time wind utvjlo wi ich are in seconds and radians respectivwlyii 


^od 

5S 

1 .0 

V a- 

*aq 

0,6 

It 

1,2 

- 0.8 

■^Itkd 

sa 

1.1 

*kkq “ 

0,8 


0.2 

^abd“ 

*ard 

as 

0.87 

*atci “ 

0.3 

*arq“ 

0.3 

^at * 

*ar 

iSS 

0,9 

^tr “ 

0,55 

acj. « 

1,1 

x^ =1.1 


sz 

0,01 

II 

0.02 

^kq * 

0.04 

= 0.001 


mam 

0.0011 

^e * 

0.05 

^e 

0,3 

^0 = ■' -0 

V 

Si 

0.019? 

Vi^ SS 

0,0032 

^0 == 

314.0 



/t tho choson operating rt’int the machine is delivering 
rated ro at unity pov«jr Tactor to bhe irjfinjto bus. Por 
the above operating conditions tho Vsicious machine 
qu'intitios aro caloulatod using the phasor diagram shown 
in Piguro 7*4 and are given below ; 

55.20 i^^« 1 ,042 3^0 » 


1,213 3.^0 
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7,4 PDnromiiiiCE am con CTTioiia HEGumTaRs 

III Lho coAivontionol design procedure, suitable 

configuraiiions for the ro£,iilators are chosen first and 

bhon bho game and oino constants arc adjusted to meet 

Lhu dosign op 00 If 1C .■'t ions and stability requirements. 

31 32 

It is cliovm * that control of quadrature axis field 
landing from load angle feedback effectively extends the 
stahLo operation to cover tlio whole range of reactive 
pov/ci"'. Por tho syston oonoidored, a voltage regulator 
\/itli one limo constant and an angle regulator witn 
single Lino cons bon I wo used for tho excitation of 
rcMOuivo iind torque v/jadings respectively. ^ conventional 
spool '^ovoaHoi* .'ith two time oonstanls is used for the 
coAAfcral or input •'or quo. In tho steady sbabo these 
eonlrollcrs will keep Uio terrain .. 1 voltage and load angle 
vil apccjJ’iod values irioupuotivo of ihc lotiding conditions . 

Tho porformance equations for bhe regulators are 
32 

given by * 


Voltage Regulators 

“ l+!t^p ^'rof *“ 
Anglo Rogulator: 



( 7 . 28 ) 


( 7 . 29 ) 
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Speed Crovcnior^: 

K 

“ ri+IgP)!T+V)^''o--''VWo (7.30) 

\/hcro 

•^r ® ^ard'^^r 

Choosing the stcilo variuitios for tho controllers as 
AE^,6Ey,A!D^ and P A the oguotions (7.28) to (7.30) 
J.VQ throv/n mbo stul/C api''co form after linearizing, as 
K„ - 

pAK^, = ,p-l« - If-AE, (7.31) 

0 . Q 

pABj . . ^ 0, AS • ^ Oji -^OgAi^ - 1 j_Ab^ 

... (7.32) 

pAl^* ASg (7.33) 

p ''''2 = ^"2 ('^• 54 ) 

Oorabinlriij; Uiuac cquabaona \/ith syotou state equation 
(7.22), the controlled s^'siiom siate model is obtained as 

X as A.J X (7.35) 

34 

'i’he various regulator pai'.mtbero ..irc solecbed as'' 

K^ = 2,0 ISr*5.0 KgsaS.O 

» 0.5 * 0.25 0?g » 0,1 

Ijj * 0.5 

The syotom matrix for tho oper'^ting point considered 

in tho e.irlicr section is given by 



0,00 1.00 0.00 C.OC C.GG C.CC C.CC C.OO C.uO O.CC C.CC 0.00 
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The rpRponso oC the oo\i,rollod system for an 
disturbance in the torque \anding field current is obtained 
hy solving' t’lp oyotciii equiuion (7.35). The response is 
shown in Figure 7.6. Tl'o oscillatory response dies domi 
in about six seconds, Evoii though the response is better 
ill an the response shown in Figure 7.5, it is still 
esc ilia Lory. An improper choice of regulator parameters 
0 . n make the system unatiiblo. 

7 ,5 CQNOLUoION 

The noud for Lbo use of on additional field 
\/iiiding to improvo thi dynamic performance of the 
synchronous fiiachinoo is briefly discussed. For d.w.r. 
synchronous naohinc, the state space model is derived 
in a linovir time Inv.Tiant fern. The machine is consi- 
dered with simple angle and voltage regulators on the 
field windings and a speed governor for c -'ll trolling the 
speed. The pcrform'aioc of Khc system with and wit]-»oui: 
convontionnl regulators io obt mod using the state space 
model. The responses arc oscillv»tory and take longer- 
times to settle down. The feasibility of using either 
opLimnl or ouboptiimj constant feedback control law 
for the d,v/,r. machine is consider od in the next chapter, 
to Improve tho systom dynamic performance* 




CT-IAPm VIII 


OPTBIAL im SUBOPT DIAI COmoi, OP 
D . .SYRCIIROFOUS GEI®R^TOR 

8.1 TU PRODUCTION 

i'll optimal otale regulator is obtained for the 
d.i/.r. oynchronous machine usin^ the state space model 
derived in the previous chapter. The optimal control 
la\/ requires tho uvailability of the entire state 
va.'iabioa for fflotisuromenbs. But in many practical 
problonio thiu io difficult i>o achieve for physical 
rotaaouo* Itcuoo a dynamic oboorver is discussed which 
i'ocoJiaLruoLu Ute utiavaila lie sbato variables from the 
av.uliiblc oui-iHit muaaun jvionlo, Tne observer introduces 
oxpoiionlialiy decaying eiTOi’ in blie reconstructed otate 
vnr1ah3 0i3 unU also tr.mssfu' functions are introduced in 
Lho fofdback paths* Mac bho estimation of state 
variables boo ones a difficult task where there are a 
largo number of state variables bub only fewer output 
vnruibleo. In such oaoos, it is most desirable to 
obtain nn optimal regulator vitoich is a function of the 
measurable outputs. If this relationship is linear and 
time invariant, then bho control can be easily imple- 
mented • 

A liiiQ ir constant output feedback control of the 
d,\/,r* synchronous maohine is also discussed. The 
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porformanco oi’ tho syatora with optimal and suboptnaal 
rof^lQl^ors and \/lLh dynamic ohserver is obtained for 
impuDno ty[)o dis^uurbanocrj* The responses are com'Dared 
nnd ooiKluoions ara drawn. 

8,2 OPTiriAL I^BGUTuVi’OR ROTi D.'/.R, FACHUTB 

Tho clas.iical morliods of control systems design 
lor s nchroiiouo n.-’.olinos assume apriori configuration 
or lihr rr'|,»iUtl.or/j. Tno feedback signal is assumed to be 
dcrivt.d tb-’oj.i ..oj.ie oiuputi quanliby; for cxaiaple, a signal 
prupoj.tio‘1 1 to iM i> ijoriiinai voluagu is fed to une field 
wiridLUtj M *inlv.iu Uic teri'iiiit'l vo' bago at desired levels. 
In jnudonn aoiiiii*ol .v^boms practioo, an inbograted form of 
cojibrol XU dv’riv.d wticii to a fu.'Ction of chc states or 
output... Tli'ia uu oi.cinol or a uuboptxiul regulator is 
ol)triined by a pi’OfWU’ ohoioo of a purformanco criterion, 

A qunUr.txc ,'urformanoe index in the suatc and 
co’jbrol V 'rxnMos lo sclootod for tho design of an optimal 
rr^'ui al.or Cor the cl.\/.r. naohinc. The problem is tuerefore 
to Cl d n CO. brol lav; wi''3ch manimizos the perfoi’manco 
index 

J « i (x'^'' a X + R u) dt (0.1) 

' 0 

subject to tho otniiu o^u Ixon 


X « A X + J1 u 


( 0 , 2 ) 
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Tli<j (jbaic model siven by equation (8.2) is derived in 
Oliapbii* VII. Tho optima] control law is shown to be”^*^ 

u - - H-’ b'"' P I (8.3) 

wi'orc P in b'ic poaib tvn definite solution of 

P /. + P - P B R“'’ r + <i = 0 (8.4) 

ind the optimally oonbaollod system is given by 

i (/' - B H'"*’ b''^ P) X (8.5) 


'I’l (• poinb eonoidored is t'lat the d.w.r. 

mueluiu t;j Ui iivt rijv,, ivtbcd KVy\. at r.jtod voibuge at unity 
powoi* f 'C'boJ* bo b fcC i*jj lube bub. For this load condition 
Lau, uybbom riubrxooo i\ and 'i arc given in the previous 
caaiJui.r. *1* lo v/v’ij^libvUf'c iin-i trices for the state and control 
V ja.iljlus a>ii uiiouon uo dioouasod in Chapter IV end are 
&vun by *4 w dia (1 ,1 , 1 ,1 » 1 ,1 ,1 ,1 ) and R » dia (1,1,1), 
i/ibh those j.abricos, the Ricoabi equation (8.4) is solved 
by the nethoU of successive approximation (diuCussed in 
Chap It r IV) and bho P-mairi:'' is obteaned as 

“1.5227 O.OOOG 0.2590 0.2054 -0.27:53 -0.0614 0.2873 0.0699 

O.OOOC 0,018B -0.0172 -0.0129 0,0184 C.0037 -0.0184 -0.0060 

0.2500 -0.0172 0.2807 0.2430 -0,3100 -0.0399 0.3249 0.0477 

U.2034 -0,0129 0.2430 0.2G16 -0.2988 0.0208 0.3136 -0,0188 

-0,?735 0,0184 -0.3100 -0.2988 0,3645 0,0114 -0,3769 -0.0174 

-0.0614 0,0037 -0.0399 0,0208 0.0114 0,0708 -0.0108 -0.0718 

0,2fJ75 -0,0184 0,3249 0.3136 -0.3769 -0,0108 0.3981 0.0172 

0.0699 -0.0060 0.0477 *0,0188 -0.0174 -0,0718 0,0172 0.0836 
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The optimally controllod system given by equation 
(8*5) IS then solved for the dynamic performance when there 
IS an iaC4:£^i disturbance of 0,05 p.u, m che torque 
winding field current. The optimal regulator response is 
shown in Figure 8.1 . The response is nonoscillatoiy and 
decays exponentially veiy fast. By comparison of this 
response with conventional regulator response shown in 
rigure 7*5f it can be concluded that the optimal regula- 
tors arc superior in performance than the conventional 
regulators . 

The optimal control law is obtained at different 
opo..*ating conditions and the average -values of the 
performance index for these conditions are calculated 
using equation (6.20)* The iiachine is delivering rated 
KVA at rated voltage to the infinite bus, at tne different 
power factors. The performance index values are given 
below for the different conditions. 


Powei 

factor 

tJ*4 IC'G 0*8 lag Unity 0.8 Icrd 0,4 lead 

4i 

J 

5.524 3.318 3.008 2.863 2.849 


From the above results it con be concluded that the 
optimal regulator obtained at one operating condition 
c--»n be used at different load conditions without much 
performance deterioration. 
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8*3 PEEPORMAITCE Op lEElE OPTBlAIi EliGDIiATCiR POH t.tttr 
EECLOSURE 

®he system shown in Pigure 7,3 is operating 
initially with one line in service and nas the following 
operating conditions in steady s-catc i 

= 65.80 Wq = 3U.0 - 0,934 i^^ = 1 .638 

^do * ^ 9,0 “ ^*^•12 v^Q = 0.76 = 1.0 

2ne problem is bo invest \"ate the performance of t-ze 
system v*ith opti lal regulc'tors, vdien the second line is 
I'ccDosed, Viti-’ two lines in service the operating 
cont'jtxons arc .Iroady given in the previous cnapter. 

Prom the control theory poinb of viev/, ■che problem is to 
ti'ansfer the system stato from one line service conditions 
to two linos operating conaitions. The performance of 
the d.v/.r, machine for this disturbance with conventional 
regulators, ducussed in the previous chapter, is shown 
in Pigure 8,2 for this sysuem state traiisfer. With optimal 
rogulabors tho performance is shown in Pigure 8,3. By 
comparing these two responsos, it can be esc*blished that 
the optimal rogulator is better than conventional 
regulators oven for such large disturbances, 

8.4 PBEPOmiilllOB OP D.W.H. MOHIIIB vTITH PYWAirEO OBBEETEF 

The optimal regulator discussed in the last section 
calls for the measurement of entire state vector for 
feedback. But it is seldom that all the state variables 
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can be measured directly. la this section, a compatible 
dynamic observer is obteinod which reconstructs the 
complete state vector from the output measurements, The 
output equation for the aj'-stem under consideration is 
derived in tlie previous chspter and is given by 

^ ° ^ ( 8 . 6 ) 
The dynamic observer is discussed in detail in Oliapter 7. 
Tor ti f d.w.r. machine, the eight state variables must be 
reconstructed from the three output variables namely A6 , 
Aw and Thus the order of tho compatible dynamic 

observer boco''s)S five. The observer matrices ? and G are 
selected such that P is a stable matrix and the pair (P,G) 
IS controllable and are given by 


10.0 0 0 0 0 
0 - 10.0 0 0 0 

0 0 - 10.0 0 0 

0 0 0 - 10.0 0 

0 0 0 0 - 10.0 


0 as 


1 

0 

0 

1 

0 


0 0 

1 0 

0 1 

0 1 

1 1 _ 



1G8 


Given the matrices iP and G as above and the s^tem 
matrices A, B and C in Chapter YII for the specified 
operating conditions, the observer matrices T, H and W 
are calculated following the procedure discussed in 
Chapter V . These are given below ; 

"0.0536 -0.0055 -0,0225 -0.0257 0.0351 0,0074 -0.0349 -0.0058. 

0.4095 0,0544 0.3992 0.2261 -0.3354 -0.0743 0.0310 -0.0193 

T= 0.3370 -0.0003 -0.0460 -0.0738 0.0155 0.1192 0.01 57 -0.0028 

0.0873 -0.0057 -0.0885 -0,0995 O.05O6 -0.1182 -0,0192 -0.0086 

5.4432 0.0541 0 0462 0.1523 -0.3198 -0.1935 0.3259 -0.2214^ 


—15. 

5700 

143.30 

-99.150 

-112.70 

44.lf|^ 



—1 . 

6380 

-9.41 

-75 .460 

-77.10 

-84.87 



__ — 0. 

2839 

2.83 

-0.014 

-29.79 

2,82J 



.000 

0.0 

0.0 

0.0 

0,0 

0.0 

0.0 

0,0 

.0 

1 .0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

.2744 

0.0 

0.0 

0,0 

0.2641 

-0.1514 

0.0 

0,0 


0,0536 -0.0055 -0.0425 -0.0257 0.0351 0.0074 -0.0349 -0.006 

0,4095 0.0544 0.3992 0.2261 -0.3354 -0.0743 0.0310 -0.019 

0.3370 -0.0003 -0.0460 -0.0738 0.0155 0.1 192 0.0157 -0,003 

0.0873 -0.0057 -0.0085 -0.0995 0.0506 -0.1182 -0.0192 -0,008 

j}.4432 0.0541 0.3462 0.1533 -0.3198 -0.1935 0.3259 -0,221_ 


Once these observer matrices are computed, tlie 
observer dynamics arc coinplotely specified. Then the 
observer is Ccucaded v/ith the optimal regulator to ooxain 


the overall control sysbcin. Ihe performance of the 
cascaded system for an disturbance of 0.05 p.u. 

xn tho torque windini,' field current is shown in Pigure- 
S .4 . The responses of the variables shown have large 
overshoot in the initial portions but decay fast 
exponentially. By choosing large negative eigen values 
for the matrix P, the error in -uhe estimated state 
variables can be reduced. But this in-croduces large 
game in the- food back paths and a compromise is 
essential belftoon thooc tv/o. It is also noted tnat the 
observer dynamics a^’o very much sensitive to the operating 
conditions and henco the overall control will bo different 
for different load conditions. Thus this type of optimal 
control will have limited applications. 

G.5 SUBOPTBiL\Ii COiJIROI (S' MOHLTE 

The opbimcil rogiiLator gives the oest performance 
for a chosen performance criterion. Unfortunately, it 
cannot be implemented in many physical problems. The 
question night, bhc'eforo, bo asked v/hether a suboptimal 
control law can bo used which is easy to implement. In 
bhis section, a control law is obtained which is a 
function of the output variables only. The control law 
is thus constrained to bo 

u a F Y 

The problem is thus modified which requires tnc solubion 
of the feedback control matrix P tnat minimizes the 
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perforriMicc indox given by equotion (8.1) subject to the 
stn'Go and output equations (8.2) uiii (8.6) respectively. 

Tile oiVoliiial feed Duck uatrix is slio 7 /n^*^ to be the 
solution oi equations (G.15) to (6.17) as discussed in 
Giiaptcr VI. Using these equations and the system 
matrices A, B and 0 for the operating point considered 
in Chapter VII, the optimal feedback matrix is obtained 
by the itci.-ativo algoritlm^*^ as 

“ C . 0011 75 0 . 002 18 - 0,002318 

-0.001728 -0.001338 -0.004625 

i"0.3'.M-1 -1,4367 1 .322 

L. 

Witli 111 1.13 t'oedb u: 1; subopt itual control matri:c f, tho - 
controlled s.yston is solved for tlie response \.lien there 
is an d i.s tur'b;nico in tho torque winding lield 

current funi tin; i.'us,;onrjc is osjovra in figure 8.5. The 
comparit3on of tliis rospOxiso \;ith the optimal response 
(figure 8.1) shov/s titiit tiie suboptimal control is good 
enough and it takes longer time than the optimal one for 
tho rospi'nisos to sot'tlc dovm and it is also nonoscilla- 
tory, liic iiiiploiauntabiun of the control lav7 is simple 
because it requires only tho output variables for 
feedback, lln.; i)orformancc deviation is not very much 
as Can bo soon by eomparison of tho two rosponses. 









CliG subopt iml control Im obtained for the above 
oporatin:. coxiclitions is used at different loads. The 
average por I’orLianoo values by use of the suboptimal 
control c.ilculatod at unity power f'^ctor and at rated 
KVA, when usod at different operating conditions are 
given below. The mchine is delivering rated KVA at 
dirforcub power xi-'ctors to the 3xifinite bus. 


iO’/Cu’ 

^ 1. ouor _ 

0.4 lag 

0#8 

Unity 

0.8 lead 

A 

J 

18.025 

17.0024 

16.96 

17.624 


Prom the wibovo results, it can be iixf erred unat blie use 
of a parbiculcT control law obtained at Oxie operating 
point can bo uood over wide rango of load conditions 
with loss porforrunce deterioration. Tho performance 
of tho ays Lon with this suboptimal control law is 
liives uigutod for largo disturbance such as line re- 
closure as di.jouGoed in Section 8.3. !I!ho performance 
obtainocl bh suboptimal control is shown in rigure 8.6. 
Comparison of this response with t?c conventional 
rogulnbor response (Figure 8.2) for ^he sene disturbance 
reve-ils Lhat suboptinal regulators are better th'n the 
conventional ones. 

^3 discussed in Chapter VI, tho measurenent of 
the output variables except the roi.or angle is easy and 
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also thexr steady state values are knovm apnori. But 
uho system rotor angle \/ith reference to the infinite 
bus changes V(ribh chc load coxiditions. Insxead of measur- 
ing 5 with respect to the infinite bus, the angle between 

terminal voltage and tho rotor quadrature sxis can be 
11 

measured . In the case of d.w.r, machine, this angle 
IS maintained constant rt 30oB by proper excitation of 
the torque and roactivo \mn(^3ngo for all operating 
cOiiditioi'S . Thus tho roi.or onglo feedback can be 
obtained ti'om tho ncaouronont of this ai^le . Then 
tho SuO. dy ot ^ to values oC A 6 and^w are zero and 

thoir x'orofonco values arc specified, fnerefore this 
subopt imal control can bo easily implencnted on 
practical systons. 

8,6 OOWOLUSION 

An optimal state regulator is obtained for the 
d.w.r, synciironous loachine using the state space model 
derived in the previous chapter. The optimality of the 
regulator at difxorcnt load conditions is then investi- 
gated, Sill CO the imploucntation of tho optimal control 
law requires tho availability of the entire stete 
variables for diroot non,itoA*ing, it is difficult in most 
oases to implement a truly optimal control. Thus tho 
ro construction of tho unav '.liable state variables by 
state reconstruction is considered. The response of the 
opt axial system with dynamic observer is obtained and 
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coEijJored bhe optimal rceponse. The observey dynamics 

ore very much sensitive bo tlie operating conditions a nd 
hence uhey are different for different load levels unlike 
the optimal control law. Also the observer introduces 
erx’or in the state variables. To obviate these difficul- 
ties, a suboptimal control is desired which is easy to 
implement* The suboptimal control law is obtained using 
uhe trace minimization procedure . The performance of 
blic systen with subopbimol regulator is then determined. 

By conporir.oa of tho responsoo, it cen be concluded that 
bhe ouboptiiiii] control is nest suitoble both from the 
poihl of inploricnbation and system dynanic performance. 

In all tho above treatment, tno problem considered 
IS a dotcrminiotic one. But in practical problems bhere 
arc always random disturbances and measuremenb errors, 
Theroforc, aiiy design should take into account the 
pres once of these random noises. In the next chapter, 
the design of optimal regulators in the presence of such 
stochastic disturbances is discussed. 



OHi'iPTEE 12 


STOOHASTIC OPTiMAIi OCiTEOL OB' 
SYECHROEOUS MOHIEB 


9,1 nTIEODUOlPIOlT 

Ine optinal control of synchronous maclaine in the 
presence of random disturbances in the state and 
uncertainty in the output measurements is discussed in 
this chap b or, a?he increase in the value of the average 
performance index is then investigated, !Phe average 
behaviour of the system in the presence of the vtiite 
noiso IS thon determined . 

fhe superiority of optimal regulators over 
conventional regulators for tlio synchronous oachine is 
well os Is-'blishod in the previous chapters, The optinal 
or subopt imal regulators required the neasurenent of the 
stabs or output, for feedback and ib is assuued that 
these variables can be measured exactly. However 
practical noasuring insbrumenbs frequently give rise to 
some errors, a'JLso the disturbaxices coning on the system 
are r^-uidon in naluro. Hence it is necess<ii:y to consider 
the offocli of those random disturbcoaces when designing 
optinal or suboptunal regulators. 
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9.2 ST. 03? STOOHaSTIO PIQBIEM 

iliQ sliOcbBslixc problem dxscussed here is blisti ib 
is roquirod to find an optimal control u as a result of 
minimization of expected value of J given by 

j (i^ I + u® E, U) at] (9.1) 

for tho linear time invariant dynamo system 

i t. A 2 + B u + w (9.2) 

y . 0 X + V (9.3) 

whoro M and v are Loro nean valued Gaussian wnite noise 
terms ixj Lno oyotom state nnd output measurements 
respect I voly. The cov nance of those noise toms ax-e 
defined as 


E(w w'^) 

w ^2 

(9.4) 

!> 


(9.5) 

m 

E(v w ) 

« 0 

( 9 . 6 ) 


^ and il-j are the wei^itage matrices on the state and 
control variables respectively. Also it is assumed tliat 
and 'Q 2 ciro positive semidefinite matrices and Ilj and 
Eg arc positive definite matrices. The noise terms are 
independent of tho initial state X(o) and the covariance 
of tho initial state is given as 

B[Z(o) *(0)’' ) = Po 


(S.7) 



TJsins bhe separation theorem^^, optimal control 
law IS derived by minimizing the performance index given 
by equation (9.1). 'I’he optimal control law is shown to 
be the s«2me as that for the deterministic system, for 
the bypo of disturbances considered in this chapter. The 
optimal control law is thus given 


u « - 1? i (9.8) 

The con broiled system and estimator states are coupled. 

The coiia'olled system is given by 

XsAX-BlX + w (9.9) 

and tJie state eubimator is 

XaAX-Bri+Sd-CX) (9,10) 

where 

I » Rlj''* B^ B (9.11) 

S = K 0® (9.12) 

B and K are the solutions of Riccati equations 

P A + a'^ B .. R^ P + * 0 (9.13) 

/ 1C + K a’^ S Rg S® + Qg * 0 


and X IS the likolyhood estimate of the state 

vector X, The moan oquaro histones of the state X and 
esbunated state X and tlieir cross correlations are 
defined as 

X = i 3 (I X?’) 


( 5 . 15 ) 
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1 


and 

X^ECXZ"-) ( 9 , 16 ) 

The solubioii of the following differenbial 
equation gives the covariance natrix J as^^ 

X = (A - B E))! + X (A - B B) + S Eg (9.1?) 

Then tio covariance of x is obtained as 

X « X + K (9.18) 

Boi sLcv'dy stcjLo avorane behaviour of the system, the 
soluulOii oT oquitAtion (S.17) co iverges bo a constant 
value, TUufj Lhi s beady obauc covc'ricincc maxriz for the 
oooiiiubov ubu'lo vector ^ is obtained by solving tne 
follov/iixs equations! v/hen the matrices A, B, B, Eg, K 
etc, are ooxislont : 

(A - B B) X + X (A - B B)^ + S Eg =: 0 (S J9) 

The increase in t ho value of the performance index due 
to the presence of noise is sho\ra to be 

AJ » tr(P (Jg + B® Eg B K) (9.20) 

Brom cqUk'bion (9*20) It is seen that the presence of 
noise (IC 0, 02^®^ increases bhe performance index 
on an average * 

9.3 SYHCHEOIIOUS M/'OHIME SYSTEM 

The system considered in Chapter II is investi- 
gated Cor liio system dynamic behaviour in bhe presence 
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of noise in the state and/or coxitrol and in the output 
measuromciits. The system maxrices A, B and 0 are given 
in Olupcor II for the operaxing conditions at 6 ** 26,3®. 
The noiso present in the state ond output are assumed to 
bo wliibe noise vath zero mean values and with constant 
covax'ianccs . The weight age matrices on the state vector 
X aijd on the control vector u arc chosen as Q-j = dia(1,1, 
1^1,1 ,],]) and * dia(1,1) resDocti-vcly. 


I’ho covtiriance maLricos Qn noise 

Lerms \i and v ai o oolocted for illustration as 
Qg = cUa(U.05 , 0,05 , 0,05 , 0.05 , 0.05 , 0.05 , 0.05) and 
Rg s aiti(0,25, 0.25, 0.25) respectively. With these 
system nnoricos blio j.iicoati equation (9.12) is solved by 
tho meuhod of sucooosive approxiiiation and the value of 
P-matrix is obtained as 


P s 


2,3747 0.0058 1.0652 -0.9802 0.9797 -0,1093 0,1198 

0.0038 0,0192 -0.0036 0.0055 -0.0047 0.0051 -0.0079 

1.0652 -0,0036 4.4356 -4.0986 4.0892 -0,0092 0.0148 

-0,9802 0.0055 -4.0986 3.7608 -3.7503 0.0104 0.0159 

0,9797 -0,0047 4.0892 -3.7503 3.7538 -0.0107 0.0161 

-0,1093 0.0051 -0.0092 0.0104 -0.0107 0.0844 -0.0811 

^0.1198 -0,0079 0.014a 0.0159 0.0161 -0.0811 0.091 1_ 


Then the equation (9.13) is solved using the same method 
for tho Steady state solution of K-matrix and is given by 



0.0573 -0.0150 

-0.1509 

-0.0150 1,1523 

-0.1041 

-0.1509 -0.1041 

0.5701 

-0.0903 -0.0746 

0.3639 

-0.0074 -0,0034 

0.0367 

0.0411 -0.0028 

-0.0960 

-0.0010 0.0262 

0.0116 


-0.0903 

-0,0074 

0,0411 

-0,0010 

-0,0746 

-0.0034 

-0.0038 

0.0262 

0.3639 

0,0367 

-0.0960 

0.0116 

0,2377 

0.0276 

-0.0556 

0.0080 

0.0276 

0.0067 

-0.0035 

0,0015 

-0.0556 

-0.0035 

0.0346 

0.0051 

0.0080 

0,0015 

0.0051 

0.0080 


Onco thu matrices P and K are obtained, the 
mal.ricos P and S can bo determined using equation (9.10) 
and (O.II) res )octiv‘ely* Then the controlled system given 
by oqu-Hiou (9 .9) IS complcbely s^jccificd. The mean square 

A 

hintoncn of the osLiraatod otabo X and their cross 
con*ola lions are oaloulated using the matrices 1 ani E. 
Then ulio sloady stato covariance mabrix x is oouained 


from equation (9*18) 

as 





'“ 0.0725 - 0.0250 

- 0,1974 

-0.1196 

-0.0111 0.4510 

-0.0025 


-0,0250 1.2140 

-0.0850 

-0.0633 

-0.0015 -0,0084 

0.0295 


-0.1974 -0.0850 

0.7174 

0.4571 

0.0481 -0,1266 

0.0155 

S= 

-0.1196 -0,0633 

0,4571 

0.2967 

0.0347 -0.0749 

0.0105 


-0.0111 -0,0015 

0.0481 

0.0347 

0.0076 -0,0058 

0.0018 


0.0510 -0.0084 

-0,1266 

-0.0749 

-0.0058 0.0411 

0.0042 

1 

-0,0025 0.0295 

0,0155 

0.0105 

0,0018 0,0042 

0.0083. 


Tho increase in tho value of performanoo index is calculated 
using oquabion (9*19) as^J =s 1 ,03122. 
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The behaviour of the optimlly controlled system 
IS then i^vcs biga Ged in uhe presence of measuremoit 
orrors only. It is assumed that there is no random 
disturbance present in the system state and/or control, 
Diffcront values of the covariances of the error vector 
V IS asQumed and the increase in the vslue of the 
performance index J are obtained as given below : 


Eg 

A J 

aia(0.1 , 0.1 , 0.1 ) 

0,1141 X 10“^ 

0.25, 0.25) 

0.3402 X 10”'*^ 

cUa(1 ,0 , 1 .0 , 1 ,0 ) 

] 

0.1732 X 10“^ 

dia(10,0, 10.0, 10,0) ! 

0.1752 X 10"^ 


The tune variation of y, , the covariance of the state 
vector X given by equation (9*18) is plotted for one set 
of values of R 2 as shown in Figure 9.1 . It is found 
that l'’c increase in the value of performance index due 
to the p-'osence of noise in the abate and/or control is 
more bhw a tJxat due to bhe measurement errors. 

9.4 COaOIiUSIOH 

I'h© oboclustio optimal control of synciironous 
machine 10 diooussod. Iho optiiaal control law remains 
the sarao with or without the presence of noise in the 
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state a.id/or output varaaUos. However, the presence 
or noise increases the porformerce inSer on an STerage. 
Henoo it, is neoessaxy to take into account such random" 
dislurhonces whilo designing optimal regulators, to 
Obtain the steady state average behaviour of the state 
variables. By separation theorem, the optimal regulator 
problem and state estimation problem are separated and 
solved independently. Iho sGcady state behaviour of t* 
sysuen is obtained in the presence of random noises. 

The error in the output measurements has lees contri- 
bution towards the inoroaee in the performanoe index 
then uho noJLso Lorm ni tho sta-cc. 



CIWPTIR Z 

ITCTHCD OTP jIMLIR'IIMt IiAJlGE DYIIAIIIC SISTEIIS 

10.1 nTi:iaDucTio?T 

' method of simplifying large dynamic system 
using SchvTura canonic j 1 form is presented. This method 
dOGS not require the computabion of eigen values and 
cigon vocuors as in 11 t,c case of other methods? thus 
rosulbi ^ HI loss computing time. The use of a reduced 
model ODoiined by utalo Vcxraable grouping technique to 
ooaioj’o] oil. OLaginal oysLorri la ^hcn investigated. The 
pcx’Lor ■'anee oj thw on.^inal and reduced models are 
corapa '•Qu » 

Lar;=;e inbtroonnuctod pouor systems demand complex 
compuu.iu lonal scliomes rcqulru^ excessive time and hence 
cost xor bhc transient and dynamic performance analyses. 
The analysos may bocono (difficult because of the limi- 
tations o'"' computer memory and Lhkj. It is usual in 
control fjjfStoins practice to analyse such high order 
sysuemo turou^ijh approximate low order systems. 

10.2 Si-MTji'jiM’ OP TIJj) PnOBIiEM 

The ijroblom posed hero is chat given a linear 
timi L.iv iri«jnt dynamic system 

X « A X + b u 
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xt IS required to find an approximate model^*^ 

X- = A* + b* u (-50,2 

such i-liou the reduced system given by above equation 
approxinately describes the behaviour of the original 
syscem given by equation (10.1). The state variables 
eliminated from equation (10.1) vail, therefore, have 
ncglif I'^le effee L on the system response; otherwise no 
simplii LC 10 uon of the sysbom would be possible. 


lOo GLiU^un'ICATlON m iiOHWAHZ CA:i01II0Ari lOTSA 

jynbum oimpiiXica oion using Schv/arz canonical 
form rcqairoo the br-inulormation of Lhc original system 
■CO Sol w.i ’-a form* The Schwarz ctuionical form is given by 

Z =s S Z + f u (10.5) 

whot .0 S IS of the form 


0 1 0 

—S'^ 0 1 

0 -S2 0 


0 0 ■ 

0 0 

0 0 


(10.4) 


0 0 
0 0 



and bhe vector f is given by 
f » (0 0 0 0 1)® 


(10.5) 
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The uraiierormation of any arbitrary form to Schwarz 
canoxjiical form i a ef fee bod by the transformation^® 

X = H Z (10.6) 

In this tlKOio u IS taken as a si.i^e variable. If u :^s 

a vector, then the transformation to Schwarz form and 
simplification by the present method will be rather 
difficult. The pair (A, b) is assumed controllable so 
that Luo Lt'-’iififoriOtUion "iven by the equation (10,6) is 
poscii'lc. TJsinc uquntjons (10.1), (10.3) and (10.6), 
the I'ol lojur.;:? roi.-itiono ax-e obbaj_icd , 

H Vi IT (10.7) 

i: f « b (10.8) 

If h^, 4i2» .♦.» hj^ represent tho column veebors of H, 
then bUo oquatious (10.7) and (10.8) can be written in 
beriAS of bhose column vectors as 


(AJi^ I Alig * -Ah^^^ : Ahjj) s (-s^b2 i ^i-Sgb^ J 


« b 


...( 10 . 9 ) 

( 10 . 10 ) 


Prom bho above two equations, tho column vectors of the 
tranoforinabion matrix H arc solved to give 

( 10 . 12 ) 

h,. Ab + 9ja 

T, , k- n-1 (10.13) 

^-k- '^^n-k^l ®n-k+1 ^n-lw+2' 



There! oro che IransTormation matruc H is completelj 
detcrmixaed if the elements s^'s are known. These are 
obtained Xi'om tiio following relations^® ; 


Bn = 


•n-i+2 

r 

n~i 


> n^-l 


■"n 


-2 

^1 


(10.14) 

(10.15) 


-Whole arc tho elements in the first column of the 

Houth n.iu-ix, dctcrininod from ohe cbaractcristic equation 
of t' t riU>m matrix A, as done for Routh stability test. 
Thus tup orniioformal-ion matrix H can bo obtained easily 
usmii Ik o>)u..l.iono (10.11) to (10.15). 


Simpl 1 . '1 cm!. ion: 

The Ox’df.r to which tho system can be simplified is 

obtained by comiuming bno elements s^j^ I'll® Schwarz 

®n-1 

matri.: 3. Tlio comparison sttirlo wiih the ratio ^ • 

n'.*2 

The coiUMu’ison is said to bo successful, if the ratio is 

greater ulian or equal bo 10. Thio figure 10 has been 

found bO 1)0 ouitabio for many probloms^^ . It depends 

upon Lc voriabloo of inb crest tc. ue retained in the 

roducud moUti aM tho ooufficionts s^*s are related to 

tho Cl ,0 ■' values of tho oystcu muorix A. This figure 

can bo doterminod by the pliysioal insight into tho 

^n“*i 

problem. Durixig bho comparison, if the ratio is 

loss thai 1, the process is tcrmi-iatedj otherwise, it is 

Op 

carriod OaI till r— is compared. If no comparison turns 
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out successful, It IS inferred txiat no simplification 
would be possible. That is the eigen values are of the 
same order of magnibudo and thus none of the variables 
can bo nc^lectod in the model. If ith comparison is the 
last successful one >10), then che system 

can be reduced to mth order where m = n-i. 


In tho reduction process, the last i-s-Cu.te 
variables arc thus assumed to navo nogligibla effect, 
if iuh comparison is the last suocosoful one. Tnen -che 
systcii equation (10,3) is partitioned into m and i 
componoxit vectors, giving 




®mm * \i 
• 


K 

m 

I ^ 

1 m 

» « • 

s 

• 


• • 

+ 1 • • 

A 

“1 


®im * ®ii 

1 

Z ! 
1 1 

i !^x 



* 

J 

, j 


u 


(10.16) 


In tho simplification of the syoucms, the last i-scste 
variables arc assumed to have less effect on the sysiem 
and thcu'oforc it is asoui*ied th'iu 





S 0 


(10.17) 


Then solving for from equation ( 10 . 16 ) 

h * “"^ii^^ajn \ 

and tlicroiore> 

4 

"m • ^mi “ii “am'“m 


(10.18) 


- + (^jn " ^11 


sr: f^)u 


(10.19) 


• • • 
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or 


S * S>^ 2* + u 


( 10 , 20 ) 


Thus the abovo equation givos tlio reduced model in the 
Sclu/arz To mi. 


To obtain the reduced model in the original fomi 
the system equation (10,1) axid bhe transformation equa- 
tion ( 10*6) arc again partitioned into m and i components 
as onown below. 



( 10 . 21 ) 


( 10 . 22 ) 


From Lho above equation, bhe voctors2^ and are solved 
in tormo of tiio voctor using equation (10.18) as 

^ "* ®ii ®im^^m * ®ii ^i (i0,23) 

and 

Zi = (iiu. - S'l S'i u (10.24) 


Thcroioro from equation (10*23)» 

^ + Hi Hmi ®li “ 


(10,25) 



where 


^ ”■ ®ii (10.26) 

Hence rroja cciuatxon (10,24) 


^ = III ^ 

(10.27) 

where 


Hi = - Hi:, Sii Si:i 

(10.28) 

and 


I?, IW HLl fl - H^i S-i fi 

(10.29) 


Prom oqu. tion (10.21), the dynamics of 2^ are determined 
as 

« 




(10.30) 

« 

Then sububiuuting for from equation (10.27), 



= (\». + V + cv 

(10.31) 

or 



i* 

« A'^ X* + b* u 

(10.32) 


Thus bha aimplifiod model using Schwarz canonical form 
is obtaa icd i.i the required form. The reduced system 
matrix ' and the input veebor are given by 


A’’ 

» 

V + Hi 

(10.33) 

and 




b’^ 

SS 

•■m \ 

(10,34) 

whore 

Hi 

, HiS; and fS arc already defined. 

/ m 1 
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10.4 Piri''Oi.‘Jl'lANCE OP THE KSIDUOED SISTELl 


To illustraxe tho metliod of simplification 
discuDSod in uho previous soctaon, the response of the 
froo ojfjuem (synchronous machine i.ifinite bus system) 
sho\vn li' i(jUre 2*1 y of Chapter 11 is considered for 
impulSG cypo distrubanoes. The free system matrix ▲ 

IS Given 1.1 Chapter II for the specified operating 
conditions \ath rotor an^e of 26, 3®. The system is 
simpliCiod by using Schvarz cdnonioal form. The forcing 
funcoLon u in oqu.itions (10.3) and (10.32) is zero. Tte 
vtcioi' b m cljofjon as b « (0, 0, 0, 0, 0, 0, 1)^ so that 
the poir (/'•, b) IS controllable. The original system 
order is Devon. It is rf<<uocd to a thud order system 
by the px'on(~nu rnc Ihod using Scb\’arz canonical form. The 
sbabc V i.‘l.'b3oa robaiiiod in tho sinplifiod model are 
A w aiil The reduced system ruarix a* is obtained 

from 0 qua 0 ion (10.33) and is givon by 


A’' 


0 

►77.72 

-0.307 


1 0 

-1.3802 -26.821 
0,5725 -0.3032 


The rosponoo of the reduced system is obtaiaed 
for dioturbanoe in the rotor angle as shown in 

Piguro 10,1, The original response of the rotor angle 
for tho onrao disturbance is also plotted for comparison, 
iProm tho u\/o responses iu is secii that the reduced model 



■jino! s,ystorr( 



by the present moxhod has negligible errors both m the 
inibial Qiii final periods of the response. Thus the over- 
all ho 'aviour of the simplified model is the same as that 
of Uio o’li.inai oysbem. 

10.5 sur.oi>i’J3La cohiiioi, asm hiducjid moebl 

file simplified mociols con be used to obtain a 

conLrol lav for the original syocem which is suboptimal. 

Jhi opcinal regulator can bu designed for the rediced model 

and ohis can bo used as a suboptimal control on the 

original system# The simplified model can be obtained 

VQXj easily for single input sysboma by using Schwarz 

canojiioal ( oi m as already discussed. Since it is 

dif’ioulu to obt-'iii a uimplifiod model by that method for 

mult IV I'j 'bio control sysboms, Lho method of state vari- 
40 

able {rouping is usod for tlio present analysis. In 
this wo I hod, uho suabo variJibles \rtuch have large time 
cons bar to nrc grouped topcbhcr as and the rest of the 
variables vhich hnvo negligible effect on the system with 
smaller wiwe constants, aro groopod together as l?hia 
typo of ,.\rouping is done i.ithor by having seme approximate 
]aio\/lcd{;o of tho oigcn values of bho system matrix or by 
physical insight into tho uacurc of tho state variables. 
Thus partitioning txic system equation 

Z * A X + B u (10.35) 

into 3^ and Xj^ component voobors and assuming that 
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« 



= 0, Lho eliminated variables are solved to give 
^1 * "^i ^im \ *“ u (10,36) 


Then ohc reduced model is obtained as 

^ ® ^Aram “ -^Di ^-11 ~ B3_)u 

or 

X s X* + B» u 


( 10 . 37 ) 

( 10 . 38 ) 


Thus rouaining AS , Aw and A i^^ js tho sbdLe variables 
Tor Uio reduced system of the oricinai system discussed 
in OliapOiT II| the reduced system matrices A'?[ and B* are 
obttn.aod vis 



0 

- 79.74 

- 0.6098 


1 

- 0,0536 

0,0003 


0 

- 22.724 

- 0.3264 


- 0.0031 0 

- 0,05 52.08 

2.126 0 


This Tiiotliod Cvin be used oven for multivariable centre! 
sysLoms '/here u can be a vector as taken in equation 
(10.35), It 10 shown^*^ that the r< sponse of the reduced 
model obtained by this method is in olcso agreement with 
the original system response in tho initial periods but 
deviates from the original in tho steady state. 3?he control 
variables used are A Tj|^ and A . 
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Subop uimal Control: 

SiiG problem posed here rs to obtain an optimal 
control lew for the roducod system given by equation 
(10#3C)jr v/hich con be used as a suboptimal control lew 
for Iho original system given by equation (10.35). The 
performance index for the reduced model is chosen os 

J = i /“ (X?*'' «* I* + u’’ K» u)db (10.39) 

Tho optiinJil control Isw for the reduced system is 
17 

obtninc d as 


u a (10.40) 

where 

P* (10.41) 

P* IS obuiuiod as n positive dofinioe solution of tho 
Riccati cqu.it ion 

P^Alf + P* - P* B* S*"’’ B*® P* + Q* = 0 (10.42) 

The weightage matrices k* and H* are chosen as identity 
matrices* Tho P* -matrix for tho reduced model given above 
IS obtained as 


P>^ « 


"l *7877 
0.0049 


Jj).10^5 


0.0049 0.1035 j 

0.0193 -0.0061 

-0.0061 0.4449 


Ifow> these throe state variables namely as , aw 
and Aj^fd aonumed to be tho output variables tnat can 
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be fodback to the origiml system. Thus the use of this 
control law^^ results in a subopimal control of the 
origt.i^a syotenu Tho response of the original 

syscoii by using ouch a suboptunal control is shoTO m 
Figuro 10.2, for field current disturbance. The response 
IS comparable with the optimal system response. It ig 
iionoscjLllabory and decays fast exponentially to the 
stcrdy state values . This control law is better than 
the convoiibional rcgulntors* 

10.6 COrOIUSlON 

Tho use of Schwarz canonicol form in simplifying 
large, linear tune inv..riant, dynamic systems is presen- 
ted, fiiLS method does not require the computation c£ 
eigen Vviluos or eigen vectors of the system matrix 
compai'od uo obhor methods and hence the tolaL 
compulxng time is much loss. If the variables of 
inocreob ore related to fast acting modes, good 
reduction would not bo possible, Tho response of the 
reduced model using Schwarz oanonical form is in closer 
agreomont with tho original system response, Henoe uhia 
mothou can be used for tho analysis of the dynamic 
behaviour of the original system by reduced models. 

Tho problem considered in the thesis is a multi- 
vuriablo system; as the method of simplificption by 
Sohwnr'j canonical form is not easily applicable for 
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such sjrstoins, the method of state varialie grouping 
1C Joed to obtain a control law which is optimal for 
tho reduced system but suboptinal for the original 
system. Tho porforaanoe of the original system is 
invootitjntod with this type of suboptimal coxitrol 
law and is found to give a reasonably good perfor- 


mance • 



CIIAP'XM 21 
COiJCLUSlOlf 


11 .1 mmAm iuro coi'jciufjxoiiS 

The voltai^G did frequsncy control of synchconous 
machines by means of opcimal regulators gives a better 
dynamic perforpauce. The control law is derived as an 
inlegrabed control which is a function of the entire 
u bilbo voct.or, Iho optimal feedback control law requires 
the neat.urowicnt of dovi aliens of the state from the post 
dLslurbnnco olcady olitc conditioiis* The deviations are 
dll a‘j cull lo ririii jn m P' icuical problems. However 
Iho nuhoptiii'.l coafcvol 'liioJ' is a function of measurable 
output vonahJea only, lo easy bo implement « 23 d can also 
bo used on the* nonlinear motiels. 

The erric'ienl use o? modern ojjtimizalion methods 
depend heavily upon the suoue variable formulation. Sor 
the synchronous 'iiachine oyslem considered in this inesis, 
a slnto model is derived m a most suitable form. Since 
the implemonlution of nonlinear controls is difficult a 
lineuriood stale space model is obtained v/ith the object 
of obtaining a iiiiour time invariant concrol law. The 
sbaoo variables ave taken as winding currents instead of 
flux linkagofj. The current variables are easier to 
measure bhun Ihe flux liniv.'i^os and also they are the 



nati’xal c^toice for interconnec bod systems ’.vberein the 
syiooOi'i p^*rf onn^noe cun oe ne^dily visualized m ternns of 
currents ra l.hcr bhun in berms of flux linkages. 

T’i-'O conventional regulator parameters have a 
marked orCect on Liic system dynamic performance aid 
thcrolore greubor leliance is placed on well designed 
control systeiiis* Existing methods (Hyquist, root locusi 
sensiwivity rnctliod, Miirovic method etc.) cf obtaining 
the best vuluoo of .,,ains and time constants in the 
coi vciioLonui regulating equipments are discussed, 3?he 
occOiid jicUiod oL‘ Xyuaunov lo uGod to obtain these 
pai'wua i. tM which ouwtu'u e- ability as /ell as better 
dyiiojiic poiM ormanco. The parametors obtained by this 
mothovl giVvT a bi btir rooponso with minimum settling time 
and Teas overahoot, 

Tho convent tonal deaif,n methods fof control 
sysi-OiUJ, na'-.uwo uprioi'i confi- urations for the regulating 
and at”ib»lizmg equipments . Tor the syotem considered 
all fell' stnbe variables irillucace the system behaviour. 
Hence an optimal control la/ is obtained wherein the 
control oignul i« a function of all the st?^tes. Such an 
optimal regulator gives a superior poi’fornanco compared 
to the oonvontionai rcgulaboro, T' e control law is 
linear and is oiwplo to implcmont provided the state 
variabJoo can bo moauiTod. ^lao the response with this 
control will be nonosoillacory. 
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The moasuromcnb of all the state variables in a 
s^sucin 10 not always posaiLlo, A compatible aynanmc 
observer bo otiiployefl bo osbinote the unavailable 
cualiO vr-i'iablas. The obaoDwer dynamos are very much 
seiisitive bo the oporain^j conditions and also exponen- 
tially decaying errors am introduced in the reconstructed 
state variables, -fliilo the optimal concrol law with 
complete state vector focdhaolc, obbaiaed for the linearized 
model about 4iti op(’riiin{ point fairly remains optimal over 
wido lortd oondj ti.ano»wibii observer cascaded and using 
OAily the output v.ii'iabloo oo reconstruct the states, the 
overall contr'ol Iviv/ iw a Xutiotion of the operating point. 
Honoo tho optinul oojitrcO law \aih all the state variables 
can be uaud evon - or ]uri;;t» dxoturbancos . The performance 
of bJio optin..liy controlJcd system for lin§ recolsure shows 
tl’at the opviiiul rurfulatoro arc superior in performance 
oomparod to couvuntiona3 rc/ulatox's. But the feedback 
system with dynjunic observer will not be optimal when used 
for such lart'c syntem d Lsturbanecs. 

Tho drawKicks of tho optimal control system cascaded 
\/ith obsci’v.r lo auoooosrully overcome by means of a sub- 
optimal control law. Tho contx'ol law is constrained to be 
a Ixioar function of the measurable output variables. The 
resulting auboptxmul control law gi'ves almost an optimal 
porformnneo* Tho pcrforn.mco deviation while using the 
particul.ff ouboptxmal control low at difforent operating 



cOiiditions la veiy loab. Thus ths suboptimai control can 
be used oven for lar^c disburbancos. This fact is estab- 
lished by invcobieatiiita the performance of the system for 
lino rocloourc'. The implementation of the suboptimai 
control 3 aw is also very simple. However it is still 
necessary to obtc^in tlr deviations of the measurable output 
variables from their post f'ult steady state conditions. 

This sL ’tc lias to bo computed before hand v/hich wall not 
be physically poosablo. So it is necessary to consider 
only thouo output vn*lil)3cu whoso steady state values are 
always known like the system frequency and machine terminal 
volta^t . Thr .mboptir'. i ooiitrol law with these output 
variviblofJ in tliun Ut' to rj lined and xhe perfomanoe is in- 
voaticsati-'d » 

The oynenrunouu aachiiio provided with an additional 
field windings h iviu/, cuLi-viuiu ouccitation systems, has 
groabur stability limiwo uvon under large reactive power 
operations* Tho optimal and suboptima i control of such 
machines is discusacdj using a linoarizod state model. 

The rotor an£;lc of juch maohlnos is maintained constant, 
liko the syolom frequency, by proper excitation of the 
field wintluifiu. Honc^c the noasuroment of the rotor angle 
for fccdb.-'olv 3S fPoicr. Thus the suboptimai control law 
which Is a Imonr function of the machine rotor angle, 
terminal voltigc and rotor speed cun bo easily implemented. 
It is eotablishod tJinl the system performance is betucr 
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W3 ch jubOiit Uial roc'alators bhan with conventional angle, 
volu'tjv^ ill* bpcjt? j'OQulalors* It is also found tnat the 
ouboi'Li^uil cuawoj law roii ins Oumc over wide r?nge of 
op^.L‘ Unj, ooacu.uic.io '.id hcnco c^nbe used even for 
lai-,"o poj'tu !’b . fci ono . 


Iho optimal cuid suboptimal controls discussed above 
asouiflc t]’S. uhe ncrouromont of erbate and output vans lies 
art' ,Koui' U'. I'o/.'V^r .'wy arror in the mcasurerients will 

I. f i cl tin " ’fija' I 'll **( 2 • Tt is established thaT presence of 
i‘ 'ltd OKI tii.it.iu*b K*t i Ja tbc ot^ibo and uncertainty in the 
lilt out mill i lat’i* ' .o / luu of the performance index on 
oil ivcr., ,v . I’l 'I i’ui'i', the dedl/;n of optimal or suboptimal 
t'oiiti'ol rfjioul <1 tile atv looount such random noises, Ihe 
pi'oliitiiii ipy liVi'Oti' 'ujon oi* otoclvistic opbimal ccmtrol of 
syiiohrtujouo i.i lOluiitiO io dlocuoecd. It lo found Tciim; the 
noioe pro .onl ui ta j^wtun abate laK^or control he a more 
contribution to tho dwtv.rioration of performance than the 
mouaurom* nt cito w. • 

Iliu cl(..»ifl;ii and unilyoia of multivariable control 
systoria Uijiiwj; th*) ot.'uo jpaco tcchniq.ucs depend largely 
on thf uoc ot dt-h-t-fc''! coiiputors ani limitations exist 
rogardinf', ooji niter lamoiv .’nd blno, Ihe fcacibility of 
nnalycing ouch ..y.-toin /ith reduced models is investigated. 
Tho ainpiaficd oynrlirorous machine model using Schirarz 
canonical lorn, /jivoo ■' un oponse \diioh is in close 
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aoreement with the original system response. !I?his method 
takes less time for 1he reduction process because the system 
eigen values and eigen vectors are not required for the 
simplification as compared to other methods. But this 
method cannot be easily applied to multivariable control 
systems. Hence a more approximate state variable grouping 
technique is used to reduce the systom size and to obtain 
an optimal regulator for the reduced sysi:em. fmis method 
IS useful even for multivariaole control systems, She 
optimal controller of reduced system txjus obtained, is used 
as a suboptimal controller for the original sysTiem. Ihe 
suboptimally controlled system response is found to be 
comparable v/ith the optimal one already obtained for 
the complete syscem. 

1 1 .2 COTOIElIBUTIOllS 

In the author’s opinion, the following are the 
contributions made to the field of syncxronous machine 
control using the modern '^nd opti.jal control theorcbio 
concepts. The listing also depicts the main theme of 
each chapter in this thesis starting from Chapter II to 
Chapter Z. 

1, An accurate stiue space model of single machine system 
IS obtained in a suitable form, 

2, A novel method of obt lining conventional regulauor 
parameters using second method of iQrapunov is presented. 
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5. An optinal regulator is obtained ^;hich is proved to be 
superior to conventional regulators, 

4, The difficulty in the measurenent of seme of the 

st ite variables is overcone by designing a compatible 
dynamic observer. 

5 , A suboptimal regulator is obtained for iiie system which 
IS oaay to implement on the nonlinear models as well. 

6, The d.w.r. synchronous macnine stale model is derived 
so that the optimal control techniques can be applied 
easily for the design of coiitrol systems. 

7, The optimal and suboptimal control of d.w.r. synchronous 
machine, which improves the dynamic performance, is 
obtained . 

8, fhe preliminary investigation of stochastic optimal 
control of synchronous n"’chines is presented. 

9, A method of simplifying linear, time invariant, dynamic 
ssrstems IS proposed end solved using Schwrrz canonical 
form. 

10, The synchronous machine control is investigated with 
a suboptimal control, obtained for the reduced model 
by state variable grouping. 

11 .3 SCOPE EOR EDROIHBR RBSE.'ilCE 

The problem investigaoed in this thesis is a single 

machine system. The analysis and design can be extended to 

multimachine systems witn no difficulties, except for the 

increase of the system size. 
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The optimal conventional regulator parameters are 
obtained by Lyapunov's *ietnod» by using a linearised system 
model. The problem can be investigated for the nonlinear 
systems. 

In the design of optmal and suboptimal regulators, 
it IS assumed that there are no constraints eitner on Ihe 
state/output or control variables. The problem can be 
studied with these constraints imposed on them. In 
particular the problem of magnitude constraint on the 
control variables resulting in bang-bang control and 
time optimal control problems can be investigaTJed. 

The problem of stocnastic control with state 
dependent noise «nd unhnovm noise statistics are worth 
studying. The stocnasxio control problem with tne 
added restricxion that the conxrol law is a function 
of the outputs only, will prove worthy of investigation. 
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APEBEDIZ A 


ROSEHBIOCK’S METHOD OF riJlTCTIOH IIIUnilZATiaiT 

The parameter optimization by Hosenbrock's hill 
climbing technique is discussed here, let p be the 
number of parameters to be optimized and 
be the variable parameters to be optimized. A set of 
p orthogonal unit vectors Hj , and a set of 

p constants e-j , 62 » * • • » are initially stored • The 
vectors 1 = 1,2 ,.,,,p are parallel to the coordinate 
axes of the p parameters. Initially the constants e^*s 
are sot equal to 0.1. TLe criterion function J is then 
evaluaied */ith k^’s set equal to zero. The value of J 
IS evaluated once more with tne parameters set equal to 
K + e^ Ti^ , where K = (k^, kg, ...» kp)^. If the second 
criterion function is smaller (success), K is replaced 
by K + e-jn^ and e^ is replaced by 5e^ s if greater 
(failure) K is left unaltered and e-j is replaced by - . 

This procedure is then repeated in turn \ath n-i replaced 
by ..., Tip, ••• until for each a 

success has been achieved and subsequently a failure. 

Then the axes are rotated in the following way. 

let the sum of all successful steps parallel to 
be denoted by d^ . The n the following sums are 
obtained : 
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«-f =5 


a« s= 


di + dg ^2 + + S '’p 

^2 "^2 + + S '’p 

• • »• ••••• 


S "p 


U.1) 


A new set of vectors is now formed by Schmidt 
procedure : 




i 

3- 

T 


= “x - I, (“3 “3 >"3 


\ “ i 


(A. 2) 


(A.3) 


for d ^ 1f2}«««yP« 


This completes one stage of "Cxie process » «7ith cnc new 
set of values as the starting, values for tx-e second sbage, 
the procedure is repeared. T.ig computa uions are 
terminated if the optimality conditioiis are satisfied. 



iI>PEl®IX B 


COBTHOIMBIDITY AND OBSERVABIIITr COlIDIiPIOIlS 
Controllability: 

A system is said to be completely state controllable 
if it IS possible to branslcr any given initial state iCt^) 
to aiy desired fin^ state Z(t^) in a finite time interval 
tjj;^ t £t^, with an tmconstrained control vector u(t), 

A system is said to be completely output control- 
lable if aic^ given initial output i(tQ) can be transferred 
to the desired output Y(t^) in finite time interval 

^ 1. » with an unconstrained con'crol vector u(t), 

Bor a linear time invariant sj^stem given by 

i = A I + B u (B,1) 

V = C X (3.2) 

the complete state conbrollability implies that "uho 

1 *5 

composite W-matriz 

V = (B : AB : A^B : ( 3 . 3 ) 

• • 

IS of rank n. 

Bor complete output controllability of uixe system 
given by equations (B.1) and (B.2) the composite T?-maori 3 : 

V (OB : CAB : CA^B : . CA^"'^B) (B.4) 


IS of rank m* 
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Observability: 

A system is said to be completely observable ua 
Che time interval t^ <. t £ t^ for any given t^ and some 
tj, if every state variable can be reconstructed from 
the knowledge of the output variables Y(t) and the 
control vector u(t) in the time interval t^ _< t ^ t^- 

Dhe system given by equations (B.1) and {B«2) 

IS completely observable if the composite '/-matrix 

w = [Q^ : : (A^)^d^ : (a^)^’''*g^] <B45) 

• • • 

IS of rank n. 

The existence of optimal control law requires 
the prior investigation of the system state and/or output 
controllability. Bor clio design of an observer system, 
the notion of observability is necessary. Tnese conditions 
are invoked before a solution to the problem is attempted. 



APPEEDIX 0 

OPIUm IIHEAE EBGUMTOE 

The linear feedback control law is derived using 
Pontryagin*s minimum principle for the system 

2 = A 2 + B u, 2(tQ) * 2 q (C.1) 

It IS required to find the control which minimizes the 
cost function 

J = ^ 2^ S 2l^^^ + ^ (2^ ^ 2 + u'”^ E u)dt (0.2) 

^ o 

ITithout loss of generality it is assumed fliat tho Q, E 
and S are symmetric matrices. The Hamiltonian is formed 
as 

H(2, u, X , t) s ^ 2^ « 2 + ^ u'*^ E u + X®(A 2 + B u) 

... (C.3) 

where \ is the Lagrangian mulciplier veccor. Application 
of the minimum principle requires that for an optimal 
control , 

^=0 = Eu + B‘^X (0.4) 

and 

3TT • 

'5J* - x=s^2 + A'^x (0.5) 

with the terminal condition that 

XCtf) . S I(-tj) (0.6) 
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Thus the optimal control is given hy 

u 5= - (C,7) 

A closed loop optimal control law is obtained by 
assuming a solution to cguation (C,5) an a form similar 
"DO the equation (0.6). Thus let 

X s= P Z (0.8) 

be a solution. Substituting equation (0.8) in equations 
(0.1) and (0.7), the following equation is obtained ? 

Z = A Z - B B^ P Z (0.9) 

and from equations (0,5) and (0,8) 

X r:PZ + PZ = -<iZ-A^PZ (0,10) 

Combining equations (0.9) and (0,10), 

(P + P A + P - ? B 3^ P + a) X = 0 (0,11) 

Since this must hold good for all nonzero Z, the term 
premultiplymg Z must oe zero. Thus the symnetric 
P— matrix must satisfy the Eiccati differential equation 

Ps-PA-A^P + PBE"'^ B^P-5 (0.12) 

\^ith terminal condition, 

P = S at t = t^ (0,13) 

Thus a closed loop optimal conbrol law is obtained as 
u = -E"'' B^ P Z (0.15) 

Tor linear time invanant systems ynth constant Q and E 
ijiabrices and for the nifiiate time regulator the 
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Eiccati differential equation (C,12) reduces to an 
algebraic equation 

PA + A% - P B E"** P + a = 0 (C.16) 

Then the control law becomes a constant feedback of the 
state vector. 
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iPIBKDIZ D 
RUFGS-m!!]?A KJIHOD 

She Eunge-Kutta fourth order integration method is 
described for the solution of -vector differential equation. 
Consider the system 

i = fU, t), z(o) = Zq (e.i) 

If the value of Z is knorm an time t, then the value of 
the solution vector can be obtained at time t + At as 

Z(t + At) = Z(t) + (K^ + 2E2 + 2K3 + (D.2) 

wliero the vectors and are calculated from 


ICj = i t f(Z,t) (D.3) 

Kg = A t f [(Z + ■^), (t + -f)Z (D.4) 

K5 = At f CU + ^), (-c + (D.5) 

K4 = At f L(Z + K^), (t + At)] (D.6) 


I he equation (D.2) is obtained by using the isruncated 
laylor series expansion of the function f(Z,t) about 
the starting point Z(t) and equating the corresponding 
coefficients of ( At)^, n = 1,2,3»4* Kie method has a 
rounding of errors proportional to ( At)^. The method 
IS numerically stable for reasonable values of -'t. 
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